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A bstract 


This thesis presents a general algorithm for the representation and approximation of 
blending surfaces between two non-uniform rational B-spline (NURBS) surfaces using 
surfaces of the same class. The algorithm permits the approximation of blending 
surfaces with up to curvature continuous boundary conditions along the intersection 
of the blending surface with the primary surfaces to be joined. These intersections 
are called the linkage curves. 

Ihe blending surfaces are approximated by low degree B-spline surfaces. created 
using a lofting method. Isoparameter curves of the approximating blending surface 
are constructed based on differential geometry properties of the underlying surfaces. 
First derivative properties are required to satisfy tangent plane continuity bound- 
ary conditions. To satisfy curvature continuity conditions. the cross-link curves are 
constructed to match the curvature of a normal section of the underlying surface at 
the end point of the cross-link curves on the linkage curves. An adaptive sampling 
procedure is used to improve the approximation. 

Finally, the algorithm is applied to the practical problem of designing a fillet for 
the intersection of a marine propeller blade with a hub. An adaptation of the method 
which creates a standard marine propeller fillet automatically is also presented. 


Thesis Supervisor: Nicholas M. Patrikalakis 
Title: Associate Professor of Ocean Engineering 
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Chapter 1 
Introduction 


The shape of objects used in engineering applications reflects both functional re- 
quirements and manufacturing processes. Assemblies may be made of parts with 
simple geometric shapes because of the manner in which they interact, such as shafts 
and bearing surfaces. or because it 1s convenient to manufacture them in that manner. 
On the other hand. complicated sculptured surfaces may be required to satisfy func- 
tional requirements such as hydrodvnamic or structural performance specifications. 
There are also aesthetic reasons for choosing a particular shape. 

When parts are assembled into the final structure. sharp edge intersections often 
arise. When these edges and intersections are internal to the assembly, or concave. 
they can cause stress concentrations or initlate flow disturbances. When they are 
external projecting features, they can pose safetv hazards. To address these problems. 
a smooth transition surface 1s usuallv provided. These secondary surfaces are called 
blends. 

In traditional design and manufacture, blending surfaces occur as the by-product 
of the limitations of the tools and processes used in their production. As such they 
are often referred to as engineering surfaces rather than design surfaces. Examples 
include the minimum radius of a cutting tool in the machining of parts. and forming 


limitations of casting, forging and pressing processes. If a smooth surface is required. 





fabrication guidelines are usually given to achieve the desired effect. Consider. for 
example, the junction of two machined surfaces. The designer can specify the size 
of the ball-end mill used to form the junction and can define check surfaces which 
bound the movement of the cutter. Á blending surface of circular cross-section will be 
produced between the check surface and the machined surface. The disadvantage of 
this technique is that the designer is unable to completely analyze the performance of 
the surface because a portion of the geometry of the surface is unavailable in explicit 
form. Further. the resulting blending surface is an interpretation of the manufactur- 
ing instructions. completely dependent on the skill and experience of the machinist. 
Modern computer aided design and manufacturing technology increasingly empha- 
sizes automation of processes and, therefore. relies on the explicit representation of 
the entire surface of an object, including secondary blending surfaces. 

This thesis examines algorithms for the representation and approximation of 
blending surfaces. A general blending surface is a surface which provides a smooth 
transition between different. but intersecting surfaces. or is one which smoothly joins 
two or more disconnected surfaces. Blending surfaces are often required on edges 
and at corners and joints. The principal characteristic of a blending surface is that it 
provides a smooth transition between the surfaces that it joins. and that it joins these 
surfaces with some specified degree of continuity. The exact shape of the interior of 
the blending surface is usually less important than the manner in which it joins the 
principal surfaces. 

While some applications simply require the presence of a blending surface. others 
may have functional requirements. such as the joint between an aircraft wing and 
the fuselage, the bulbous bow and hull of a ship, or the blade and hub of a marine 
propeller. When the blending surface plays a part in the ultimate performance of 
the component. it becomes a design surface which must be accurately specified in the 
design process. If the surface is left as a traditional by-product type of surface. the 


designer can not guarantee the performance of his design. 
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Having established the need for the explicit representation of blending surfaces. 


the general characteristics of such a method are considered. The method should: 


e be applicable to a wide range of geometric shapes and configurations. It must 
be able to join simple planar surfaces. algebraic surfaces such as quadrics. tori. 


and surfaces of revolution. and complex free-form surfaces. 


e be predictable in its use so that it can be automated. Given that. it must 
also offer sufficient flexibility so that the designer can manipulate the blending 


surface in regions where problems arise. 
e allow tangent plane continuity along a common boundary between surfaces. 
e allow surface curvature continuity, 1f required. across a common boundary. 


e allow comparison of the surface position and differential properties along the 


common boundary for verifying continuity of the required degree. 


e generate blends of a local nature. The method should not alter the global 


geometry of the surfaces to be joined. 


e allow high accuracy approximation to make possible the integration of blend- 
ing surfaces in modern geometric modellers and to address the need for high 


accuracy manufacturing operations. 


[n this thesis. a general method for the representation and approximation of blend- 
ing surfaces which join two non-uniform rational B-spline (NURBS) surfaces with a 
surface of the same class is implemented. NURBS surfaces were chosen to repre- 
sent the primary and the blending surfaces because the form provides a general. 
unambiguous mathematical representation able to represent a wide range of surfaces. 
Additionallv, there are complete standard data formats for NURBS surfaces which 
provide for the direct exchange of geometrical information between different geometric 


modelling, design and manufacturing systems [Smith 88]. 
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To better understand the method. a brief description of some useful terms 1s 


provided: 


underlying surfaces. Two or more surfaces which are to be joined in a smooth 


manner. These surfaces are also known as primary surfaces. 


linkage curve. A curve which defines the extent of the blending surface on the 


surface to be joined. The linkage curve is the intersection of a primary surface 


and the blending surface. 


cross-link curve. A planar curve which connects corresponding parameter values 
on each linkage curve. The cross-link curves form a family of curves from which 


the blending surface is lofted. 


directional curve. À curve in three-dimensional space which determines the orien- 


tation of the cross-link curves, given the endpoints determined by the linkage 


curves. 


bias functions. A control parameter which the designer may use to influence the 
shape of the interior of the blend. Bias functions typically affect the first partial 


derivatives of the blending surface across the linkage curves. 


boundary conditions. The degree of continuity with which the cross-link curves 


and the resulting blending surface join the underlying surfaces. 


Figure 1.1 illustrates the relationships of the supporting curves. 

The algorithm approximates the blending surface with an integral B-spline surface 
which is up to curvature continuous to the underlying surfaces. The linkage curves are 
approximated to within a specified tolerance and the knot vectors are merged. The 
cross-link curves are created using the common knot vector of the linkage curves. the 
directional curve, and the parametric partial derivatives of the underlying surfaces. 


up to second order, if required by the boundary conditions. The cross-link curves are 








CURVES 


CURVE CROSS-LINK CURVES 


Figure 1.1: Components of a blending surface 


then lofted to form the blending surface. An adaptive sampling procedure is used to 
test the blending surface against specified tolerances. Cross-link curves are added, 
where necessary, and the surface relofted until the boundary conditions are satisfied. 

A significant contribution of the algorithm is the degree of control provided to the 
designer over the boundary conditions along the linkage curves. Most blending algo- 
rithms permit position and tangent plane continuity across the boundaries. Position 
continuity 1s important for ensuring the blending surface in fact joins the underly- 
ing surfaces. By itself, however, position continuity cannot ensure smooth transition 
surfaces. Tangent plane continuity is the most common type of boundary condition 
requirement, and is adequate for most applications. Indeed, [Mortenson 85] observes 
that most mechanical parts have rounded edges which blend into planar faces for 
which tangent continuity is sufficient. For more complicated surfaces such as sculp- 
tured or free-form surfaces, there can exist significant curvature discontinuities along 
such a boundary. For aerodynamic and hydrodynamic surfaces, curvature discon- 


tinuities can lead to undesirable discontinuites of fluid particle accelerations. The 
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algorithm given in this thesis allows the accurate approximation of curvature con- 
tinuous blending surfaces for both simple and complex underlying surtaces. One of 
the features of the algorithm we developed. contributing to its flexibility. is that the 
boundary condition along one linkage curve is independent of the boundary condition 
specified along the other linkage curve. It is possible. for example. to represent a 
blending surface which is curvature continuous along one boundary while being only 
position continuous along the other. 

The creation of complete geometric models of complex structures encountered in 
marine applications has been an area of particular interest in the MIT Ocean En- 
gineering Design Laboratory. For an application of the algorithm we examined the 
creation of blending surfaces between marine propeller blades and hubs. Increasingly 
demanding performance requirements coupled with more complex shapes are forc- 
ing marine propeller manufacturers to use numerically controlled (NC) machining in 
place of traditional hand-finishing methods. While the use of NC machining has the 
advantage of repeatability between components, improved performance. and reduced 
unit production costs, it requires an explicit mathematical representation of the en- 
tire component. An area of particular concern in the design and representation of 
marine propellers is the junction of the propeller blade with the hub. This sharp 
edge intersection presents both structural and hydrodynamic concerns which are ad- 
dressed with a feature called a fillet. The fillet is simply a blending surface which 
contributes to the structural strength of the blade and helps to reduce vibrational 
blade loading, while providing a smooth transition surface which joins the principal 
hydrodynamic surfaces to reduce adverse hydrodynamic effects. Several examples of 
marine propeller fillet blends are given to demonstrate our algorithm. 

The thesis is organized as follows. 

Chapter 2 presents a survey of previous work in the area of blending surface 
representation and approximation. 


Chapter 3 presents a formulation of the blending surface problem addressed in this 
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work. The presentation begins with an overview of the basic steps in our algorithm 
for the design and representation of a blending surface. Various types of surface repre- 
sentation are considered. and our reasons for the use of non-uniform rational B-spline 
(NURBS) surfaces are given in some detail. A brief review of integral and rational 
B-spline curve and surface theory ıs included. The chapter continues with details 
concerning the specification of linkage curves. directional curves and boundary con- 
ditions. As noted above, the algorithm allows boundary conditions along the linkage 
curves of position, tangent plane and curvature continuity. A discussion of certain dif- 
ferential geometry properties of curves and surfaces necessary for the representation 
of curvature continuous blending surfaces concludes the chapter. 

Chapter 4 describes in detail the implementation of our method. including the 
approximation of the linkage curves in terms of low-degree B-spline curves, the con- 
struction of the cross-link curves, and the approximation of the blending surface by 
lofting of the cross-link curves. Improvement of the blending surface through adap- 
tive sampling for tangent plane and curvature continuous boundary conditions is 
discussed. 

The method is illustrated with a series of examples in Chapter 5. The examples 
demonstrate the importance of proper parameterization of the linkage curves. Control 
of the shape of the blending surface with the bias functions is illustrated. along with 
two biasing methods. We end with several examples of marine propeller blends. which 
serve to illustrate the practical applicability of the method in marine propeller design 
and manufacturing. 

The thesis concludes with a summary of the method and recommendations for 
further work in Chapter 6. 

Two appendices are also included. 

Appendix A provides general information on marine propeller geometry and defi- 
nitions of the relevant terminology to aid in understanding some of the examples. 


Appendix B gives the details for our implementation of an automatic blending 





technique for the standard 3T-T/3 marine propeller blade-to-hub filet. 





Chapter 2 


Survey of Previous Work 


2.1 Introduction 


Various techniques for generating blending surfaces have been reported in the liter- 
ature. [Woodwark87] presents a thorough survey of earlier methods, the simplest of 
which involve superficial blends. These methods consist mainly of adding attributes 
to the description of an intersection between two surfaces. Since there 1s no explicit 
representation of the blending surface. the methods have limited design value and. 
in general, are only useful in providing simple manufacturing instructions. Early 
techniques which do provide explicit representation of the blend surface are polyhe- 
dral blending methods. Smooth edges are created on polvhedrally-defined objects 
by using recursive subdivision methods to chop-off corners. or by using weighting 
factors to make local topological changes. These are repatching techniques where 
parts of a surface are replaced with smaller patches. Surface blending techniques ex- 
tend these methods and involve the generation of free-form surfaces to smoothlv join 
other surfaces. The original surfaces may or may not intersect. The techniques are 
primarily implicit formulations based either on Liming’s method of cross-sectional 
design using conic sections {Faux 79], or on potential methods such as that given 


by {Hoffmann 87c]. Volumetric blending techniques are common to Constructive 
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Solid Geometry (CSG) systems and Boundary-Representation (B-Rep) solid mod- 
ellers |Rossignac 84] [Rossignac 85]. These methods involve the creation of new sur- 
face geometry and restructure the model to incorporate topological properties to 
permit further set operations on the object. The models may be either parametric or 
implicit. 

The following sections will discuss in some detail several current methods proposed 
for blending surface generation. The discussion begins with spherical and circular 
blending techniques. including the rolling ball blend. a popular technique for describ- 
ing blending surfaces generated during machining (Section 2.2). Parametric surface 
blending techniques which more easily deal with sculptured and free-form suríaces ap- 
pear in Section 2.3. The techniques are further divided into tangent plane continuity 
methods (Section 2.3.1) and methods offering higher order continuity (Section 2.3.2). 
Recent work in representing blending surfaces as solutions to boundary value prob- 
lems posed in terms of partial differential equations is reviewed (Section 2.3.3), as 
well as energy-based variational methods of free-form surface design applicable to 
the blending problem (Section 2.3.4). Finally, as a particular example of parametric 
surface blending, some information is presented on the design of blade-to-hub fillets 


for rational B-spline (NURBS) representations of marine propellers (Section 2.4). 


2.2 Spherical and Circular Blends 


As noted in Chapter 1. blending surfaces have traditionally been the by-product of 
manufacture, not design. Several techniques appear in the literature to model blends 
which would be produced by ball end mills or disk cutters on numerically controlled 
(NC) machines. If one can model the blend in this way one might then hope to design 
and represent it in a more explicit fashion. 

|Pegna 87| and [Pegna 88| present methods for spherical and circular blends based 


on sweep theory. Kinematically, the blend might be represented as the envelope ot 





a sphere moving while in contact with the underlying surfaces. [he analogy is then 
extended to sweeping a circular arc. The method is an application of generalized 
cylinders. and is derived from rigid sweep theory, a class of three-dimensional variable 
sweep geometric models (Pegna 87]. À swept surface is created by a planar curve or 
cross section. fixed in the normal plane of the Frenet frame of a curve (Farin 88], which 
is then carried along a twice differentiable space curve without stationary points or 
points of zero curvature. The cross section curve is called the generatrix and the 
curve which defines the frame’s behavior in space is the directrix. Reduced rigid 
sweeps are derived from this basic swept surface model by relaxing the requirement 
that the generatrix be planar. The frame's motion is contrained to be a type of rigid 
body motion known as Darboux's motion. a class of motion which occurs naturally 
in the study of systems of orthogonal curvilinear coordinates. Darboux's motions 
constrains the rotation of the Frenet frame as the origin of the frame translates along 
the directrix. À more general class, defined as rigid sweeps, decouples the rotation of 
the frame from its position, thus relaxing the Darboux rigid body motion requirement. 
This technique defines the motion of the frame with two directrices, one to describe 
the translation of the origin of the frame, and the other to give the orientation of 
the frame relative to the translational trajectorv. To complete the description. rigid 
sweeps may have a variable generatrix. one which deforms as the frame moves. 
[Pegna 88] proposes several blending techniques using rigid sweeps. These tech- 
niques are referred to as either spherical or circular blends with constant or variable 
radius. Constant radius spherical blending describes the envelope of a sphere with 
constant radius E moving in contact with the surfaces. The path traced by the center 
of the sphere is the directrix and is the intersection of the normal offset surfaces of 
the two underlying surfaces at the same offset distance. The generatrix is an arc of 
a great circle of the sphere. The same surface can be generated using a constant 
radius circular blend. For a circular blend, representing a disk cutter, the orientation 


of the disk must be controlled to maintain tangency with the underlying surfaces. 





This additional control represents the difference in a blend machined with a ball end 
mill on a three axis machine and that of a disk cutter on a four or five axis machine. 
Additional blending methods proposed allow the radius of the sphere or disk to vary. 
The methods require the directrix to be C! continuous and that it be equidistant from 
both underlving surfaces. For variable radius spherical blending this distance 1s the 
radius of the sphere. [Pegna 87] provides detailed results of these blending techniques 
and closed form solutions where possible. The blends produced ensure continuity of 
position and tangent plane along the linkage curves. The implementation is proce- 
dural since 1t requires the evaluation of the intersection of the offset surfaces to find 
the directrix, and the orthogonal projection of points onto the underlying surfaces to 
find the boundaries of the generatrix. 

[Sanglikar 90] present a mathematical model for the rolling ball blend. The 
blended surface as well as the underlying surfaces are parametric surfaces. The ap- 
proach is called modulation and is similar to the spherical and circular sweep blends 
of [Pegna 87] and [Pegna 88]. The technique splits the blending problem into the 
design of two spatial curves, a spine and a cross section or pattern curve which corre- 
spond to the directrix and generatrix, respectively. To form the modulated surface. 
the pattern curve is dragged or swept along the spine. The pattern curve may or 
may not deform as it moves. If the pattern curve is planar, the orientation relative 
to the tangent of the spine may be varied to satisfy different boundary conditions 
or other user-specified relationships. As with Pegna’s methods. the surfaces to be 
blended are offset a distance R along the normal vectors to the surfaces. with each 
surface offset the same distance. The intersection of these surfaces forms the spine 
of the blend surface, and traces the path of the center of the rolling ball of radius R. 
Since the offset surfaces are defined in terms of the underlying surface parameters. 
the intersection curve can also be expressed in terms of these parameters. 

To define the pattern curve, a planar curve is generated which represents an arc 


which is a portion of a great circle of the rolling sphere tangent at its endpoints to 





the underlying surface. The set of resulting points of tangency represent the linkage 
curves referred to in other blending methods. Linkage curves. however. are never 
explicitly defined. Rather. the method generates a family of planar pattern curves. 
each of which is circular. of radius R. and of variable arc length. A mathematical 
expression is then derived for the blend surface. expressing the patch in terms of the 
spine curve parameter. the normal vectors to the surface at the linkage point expressed 
in terms of the spine parameter, and a shape factor h. 0 « h < R. For sculptured 
surfaces. however. the resultant expression may be of high degree. although not as 
high as an equivalent implicit model. As with [Pegna 87] and [Pegna 88]. the method 


is procedural because it requires the evaluation of an offset surface intersection. 


2.3 Parametric Surface Blending 


2.3.1 Tangent Plane Continuity Blending Methods 


|Varady 89] discusses a technique for blending parametric surfaces which has been im- 
plemented in the solid modeling system FFSolid. The system operates on Boundary- 
Representation (B-Rep) models and uses quadratic curves and double quadratic sur- 
faces. The blending technique is implemented as one of a set of operations that 
may be performed on a model and generates an explicit blend surface. Ihe blending 
surface may be operated on by any of the other solid model operations. including 
reblending operations. 

The technique involves the replacement of an edge or a sequence of edges by a face 
or sequence of faces. This creates a smooth surface which joins the faces on either side 
of the replaced edge. The set of operations includes the creation of linkage curves on 
the faces adjacent to the edge, and then trimming the surfaces to the linkage curve. 
The blending surface is generated by fitting a surface joining the underlying surfaces 
at the linkage curves with position and tangent continuitv. When a sequence of edges 


is being blended, the technique permits sharp edges such as those which arise at a 
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corner to remain or the edge may be smoothed. To create linkage curves a series of 
points on each surface are specified. Án adaptive quadratic curve fitting technique is 
then used. The series of points may be selected to reproduce many tvpes of blends. 
including pseudo-rolling ball blends. In general. a point on the edge is selected and 
the plane normal to the edge at that point is constructed. Points on the surfaces can 
then be identified which also lie in this plane and are at a specified distance from the 
point on the edge. The distance can be constant or a function of the position of the 
point along the edge, and can be measured either as arc length along the intersection 
curve of the plane and surface. or as an offset distance from the point on the edge. Á 
sorted data base ıs created which provides a position and corresponding unit tangent 
vector for every point relative to a base datum. The linkage curve is then fit to the 
data points using a least square quadratic curve segment fitting algorithm. 

To complete the blend, quadratic curve segments are generated between corre- 
sponding points on the linkage curves using either a normal plane or an arc-length 
correspondence. This set of segments makes up a double quadratic strip, a lengthwise 
double quadratic surface one patch wide, and explicitly defines the blend. The blend- 
ing surface is defined using the same type of geometric elements used by the other 
operations in the B-Rep modeller. All solid model operations may be performed on 
the blend surface. including reblending. 

[Bardis 39] presents a technique for constructing tangent continuous blends be- 
tween non-uniform rational B-spline (NURBS) surface patches. using either open or 
periodic linkage curves. The method approximates the blending surface by a non- 
uniform. bicubic integral B-spline patch which interpolates position and normal vector 
on the linkage curves at an equal number of points on both patches. They use two 
methods to construct linkage curves. The first method constructs rational B-splines 
in each patch parameter space by offsetting the intersection of the patches on each 
patch using a variable distance along the geodesics (see section 3.3). 


Another method involves specification of the linkage curves arbitrarily bv the 
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user as proposed by [Hansmann 87]. A further variant implemented is a combination 
of the these two techniques. Once defined. the linkage curves are approximated in 
three- dimensional space with a common knot vector. Using this same knot vector. 
the unit normal vectors to each linkage curve. which are at the same time tangent 
to the respective surface, are approximated bv integral cubic B-spline curves. The 
approximation technique is an iterative process which ensures convergence to within 
user-specifled tolerances. The resulting blend surface matches position and tangent 
plane continuity along the linkage curves. The shape of the blend surface may be 
controlled with user-specified bias functions. The bias functions are applied to the 
unit normal vectors to the linkage curves. which also lie on the tangent plane of the 
primary surfaces. and thus control the position of the control polvgon vertices. 

[Koparkar 91] presents a technique of curve modulation for the blending of para- 
metric surfaces. Curve modulation is similar to rigid sweep theory. The important 
features of general modulated surfaces are discussed with the citation for rolling ball 
blends by [Sanglikar 90]. [Koparkar 91] modifies the modulation technique to design 
an appropriate pattern curve at each selected point as one traces along the directrix. 
The pattern curve is designed to satisfy position and tangent plane continuity re- 
quirements on the underlying surfaces. The tangent directions are selected as either 
perpendicular to the linkage curve. or tangent to the surface in a plane established 
by corresponding points on the two linkage and the directrix. The parameterization 
selected for the pattern curve is a rational form of a conic arc. A conic shape factor is 
used to control the blending surface as either elliptic, parabolic. or hyperbolic. The 
shape factor can further be made a function of the directrix parameter w such that 
the blend varies along the directrix. 

Two techniques are proposed for the design of the directrix and the linkage curves. 
In the first method. the surfaces to be blended are offset along the surface normal 
vectors a distance A, and Ay. The intersections of the offset surfaces form the di- 


rectrix. This is similar to the spherical and circular blending methods [Sanglikar 90]. 





[Pegna 37], and [Pegna 88] but the offset distances R, and R2 are allowed to differ. 
In the special case where R4 and R; are equal and the pattern curve is circular. a 
spherical or circular blend in achieved. The linkage curve is derived automatically 
using a geometric correspondence between the normal offset surfaces and the under- 
lving surfaces. À point (u. v) on an underlying surface is mapped to the same (u. v) 
parameters on the normal offset surface. The intersection of the two offset surfaces is 
a curve which can therefore be expressed in terms of the parameter values of either 
surface. The linkage curve 1s defined as the mapping of the intersection curve onto 
the underlying surface using the correspondence of the common parameters. 

The automatic generation of the directrix and linkage curves from the user-defined 
offsets at distances R, and R$ may not be sufficient for more complicated applications. 
To provide the designer with control over the shape and position of the linkage curves. 


a method employing fanout surfaces is used. Unlike the offset surface. which depends 





only on the original surface, the fanout surface depends on both the surface and the 
linkage curve on that surface. The fanout surface is constructed as a ruled surface 


described by 
D^ sDESPu(s),v(s)-t- P(u(s), v(s)) 


where P" is a point on the fanout surface with the parameters s and £. P is a point 
on the underlying surface with parameters u(s),v(s) describing the linkage curve in 
the parameter space of the surface. and P is the unit normal vector to the surface 
at P. By holding the parameter value s constant and varving t. isoparameter lines 
along the surface normal vectors are obtained. The other set of isoparameter lines 
is formed by holding t constant and varying s. Two fanout surfaces are constructed 
from the user-specified linkage curves from both surfaces. Once again. the spine is 
the intersection of these surfaces. The linkage curves are the orthogonal projections 
of the directrix onto the underlving surfaces. 

The methods presented by [Koparkar 91] are procedural. The offsetting technique 


works well for simple blends where automatic linkage curve generation is acceptable. 
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For more complex geometry the fanout surface method permits greater designer con- 
trol but at the expense of greater computational complexity. 

[Tuohy 91b} and [Tuohy 91c] present a method for the approximation of a blending 
surface as an integral bicubic B-spline surface which is tangent plane continuous with 
the underlying surfaces along the linkage curves. The method requires linkage curves 
specified on both surfaces and a third space curve called a directional curve used to 
control the orientation of the construction plane of the cross-link curves. The cross- 
link curve establishes an isoparametric correspondence between the linkage curves. 
and is similar to the generatrix or pattern curve in other cited methods. [he linkage 
curves are first approximated with integral cubic B-splines curves. The cross-link 
curves are then constructed such that they are position and tangent plane continuous 
to the underlying surfaces. The cross-link curve also hes in a plane containing a 
point on each linkage curve and a point on the directional curve. The magnitude 
of the tangent vectors at the ends of the cross-link curve is controlled with bias 
functions. The bias functions control the shape of the cross-link curve and can be 
defined independently for each boundary. Once the cross-link curves are defined. they 


are lofted to create a rational B-spline surface. The process is iterative and converges 


to a user-specified tolerance. 


2.3.2 Curvature Continuous Blending Methods 


[Hansmann 85] and [Hansmann 87] present the fundamentals of a general. interactive 
method for the design of smooth transitions or free-form blending surfaces between 
general surfaces with curvature continuity. As examples of free-form surfaces. inte- 
gral cubic tensor product B-spline patches are used. For analytically defined surfaces. 
quadric surfaces parameterized in terms of transcendental functions are used. Linkage 
curves are cubic B-splines interactively defined in the parameter space of the surfaces 
to be joined. Differential geometry is used to develop tangent plane and normal cur- 


vature properties of the primary surfaces along the linkage curves and in a direction 
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orthogonal to these curves. These properties. along with the linkage curves. provide 
the boundary conditions for the connecting surface which is a lofted parametric sur- 
face using cubic B-spline blending functions. The technique permits up to curvature 
continuous transitions across the blending surface boundary. The blending surface 
must be treated as a procedural surface. however. since the surface is irrational with 
radicals of high degree polynomials. 

[Pegna 89a] presents a criterion for guaranteeing second order smoothness or cur- 
vature continuity for blending surfaces. Curvature continuity across a linkage curve 
generally means that the normal curvatures agree in every direction at every point 
along the linkage curve. There are many ways of formulating this situation. and some 
of these are discussed in Chapter 3. The contribution in the paper by Pegna and 
Wolter is the Linkage Curve Theorem. This theorem states that two surfaces joined 
with first order smoothness, or normal continuity, along a first order continuous link- 
age curve can be shown to be second order smooth if the normal curvatures on each 
surface agree in one direction other than the tangent direction to the linkage curve. 
There is no requirement on the direction other than it not be tangent to the link- 
age curve. Furthermore, since curvature continuity is achieved by matching normal 
curvatures in just one direction, the problem of constructing second order smooth 
blending surfaces is effectively reduced in dimension to one of curve fairing. 

Using these results, [Pegna 39] presents a method for designing blending surfaces 
which ensure curvature continuity. The method is an extension of Pegna's work with 
rigid sweeps cited earlier [Pegna 87] [Pegna 88]. Here, a space curve is introduced as 
the directrix and is the axis of a rigid sweep. The linkage curves are constructed as 
orthogonal projections of the directrix onto the underlying surfaces. An isoparametric 
mapping is used to relate the two linkage curves and the directrix. Because the linkage 
curves are derived from the directrix the problem of undesirable twist in the blend 
caused by poor parameterization of independently designed linkage curves is reduced. 


To create the blend a swept surface is used with a generatrix which has endpoints on 





the linkage curves, end tangents parallel to the surface tangents in the plane formed by 
the isoparameter mapping, and end curvature such that the normal curvature of the 
underlving surface and the blend is identical at the endpoint. again in the direction 
formed by the isoparameter mapping. As examples. a generatrix of Bézier form is 
used. with third. fourth and fifth order solutions demonstrated. The fourth and fifth 
order solutions have design parameters to control the steepness of the blend at the 
linkage curves. [Pegna 89] shows that the blend is curvature continuous if the directrix 
and the underlying surfaces are second order continuous and the design parameters 
are first order continuous functions of the directrix parameterization. Furthermore. 
the isoparametric mapping ensures the end tangent vectors to the generatrix are 
not tangent to the direction of the linkage curve, thus satisfving the Linkage Curve 


Theorem of [Pegna 89a] to ensure curvature continuity. 


2.3.3 Partial Differential Equation Blending Methods 


A method for the representation of surfaces as the solution to partial differential 
equations and its application to the creation of blending surfaces are described in 
[Bloor 89a] and [Bloor 89b]. The blend is defined as a function which solves a bound- 
ary value problem. The boundary conditions are typically the position of the blending 
surface on the underlying surfaces (the linkage curves), and a number of derivatives on 
the boundary as determined by the degree of continuity required between the blend 
and the underlying surfaces. The blending surface is represented by the parametric 


function X = X(u.v) where X is the solution to the partial differential equation 


LM (X) = F(u,v) 
where L;?, is a partial differential operator of order .M in the independent parameter 


values u,v. The vector valued function F is generally zero for blend surfaces. The 


class of operators used are elliptic partial differential equations (PDE). The equations. 
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an example of which is the minimal surface equation. are well-behaved. smoothing 
functions which generally possess unique solutions without local oscillations. For 
blending surface design this leads to smooth surfaces without self-intersections or 
inappropriate intersections with the underlying surfaces. 

In [Bloor 89a], a second order PDE is used to create a position continuous blend. 
Io control the shape of the blend between the linkage curves a design parameter is 
added to the PDE. The parameter permits varying the influence or the smoothing 
effect of the operators in the two parametric directions. lo achieve tangent plane 
continuity a fourth-order PDE is given, also with a design parameter for control. 
Higher order continuity 1s suggested (a sixth-order PDE for curvature continuity} 
but not demonstrated. The examples presented involve analytic solutions involving 
transcendental functions and as such have limited use in sculptured surface design. 
Numerical solutions are also examined and these are expected to have a greater range 
of applicability. 

In [Bloor 89b], the role of the design parameters in the PDEs is examined more 
closely. The nature of the elliptic PDE allows derivation of closed form solutions 
to the equations provided the boundaries are properly selected. Design parameters 
included in the solution permit the designer to control the magnitude and direction 
of the tangent vector at the boundary. A smoothing parameter controls the relative 
smoothing of the dependent variables between the parametric directions. The effect 
is to control the extent of propagation of boundary effects into the interior of the 
surface. Both design parameters can be either constants or functions of the surface 
parameters. As functions they can be specified at either or both boundaries. and can 
be either in phase or out of phase. 

In [Bloor 89c], the method is extended to the generation of n-sided patches. 
| Woodwark87] briefly discusses topological problems associated with sculptured sur- 
face engineering. Many common surfaces cannot be represented bv a regular array of 


rectangular patches. A polvgonal hole will normally appear somewhere in the surface. 





Techniques which subdivide patches to solve this problem tend to relocate the hole. 
but it still appears. The hole generally has n sides. where n is an integer not equal 
to 4. As before. the blending surface is the solution to a boundary value problem with 
elliptic PDEs. The significant adaptation is the parameterization of an n-sided region 
with two parameters. The boundary conditions are carefully selected to include more 
than one edge in Euclidean three-dimensional space into the isoparameter edges of the 
two-dimensional mapping. Junctions are represented by singular points with certain 
characteristics depending on the behavior of the underlving surfaces and their inter- 
sections. The set of equations requires a numerical solution. The blending surface 
which results joins multiple surfaces and edges with a single parametric patch with 
no self-intersections, and. depending on the order of the PDE. satisfies any required 
degree of continuity at the boundaries. 

[Bloor 90] shows the PDE method may also be used for the design of free-form 
surfaces. In this case the boundary conditions are used to manipulate the surface 
since the designer is free from boundary and continuity conditions imposed on him 
by the blending surface problem. Internal surface features which cannot be created 
from the boundary data directly can be introduced by subdivision of the surface. 
introducing a new boundary to provide the desired feature. An example is the hard 
chine of a ship hull. A successful implementation of the PDE method would require an 
interface which translates a less than intuitive PDE design parameter manipulation 
into a natural design activity. In the case of blending surface design, the method may 


be useful for the automatic generation of blends. 


2.3.4 Energy-Based Variational Blending Methods 


[Celniker 89] proposes a method of free-form surface design using energy- based surface 
models. The method takes advantage of the physical fact that objects deform. when 
subjected to loads and constraints. to a shape which minimizes the stored energy of the 


deformation. In this respect. the method is similar to PDE methods. Surface shape 
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is controlled bv boundary conditions. external loads distributed on the surface. and 
an energy functional. Functionals are integrals which assume maximum or minimum 
values and are the basis of the calculus of variations [Fox 387]. For energy-based 
surface models. an energy functional is specified which depends on surface shape. 
[Celniker 89] implements the technique with parametric surfaces of the form x(a) = 
[x1(a), z2(a), z3(a)], where a = [u. v]. The mapping is allowed to change with time. t. 
such that the shape of the surface is described by x = x(a,t). The energy functional 
selected in the work depends on the shape of the surface and vields a scalar value for 
the energy of deformation. 

Energy-based methods seek to minimize the internal energy of a svstem. Bv 
defining a functional based on surface curvature and elongation. surfaces deform to 
naturally fair shapes. Boundary conditions and position constraints may be applied 
to any point on the surface. The choice and number of constraints applied at a 
point give the user the means to define precisely the shape of the surface in both a 
local and global manner. The shape of regions with free points (points with no or 
few constraints) is controlled bv application of forces or loads. The implementation 
uses forces such as gravity. springs, and pressure as physically meaningful forcing 
functions. Surface fairness is maintained under external loading by the internal forces 
of the surface. Internal forces are found bv taking the first variation of the energy 
functional and setting it equal to zero. 

The control of shape by minimization of surface energy might be applied to the 
design of blending surfaces by constructing a parametric surface between linkage 
curves on the surfaces to be joined. and applving appropriate boundary conditions. 


The shape of the interior of the blend could be easilv controlled bv pressure loading. 





2.4 Marine Propeller Fillets 


Marine propellers represent a particularly complex application of sculptured surface 
design. The junction of the propeller blade with the hub presents both structural 
and hydrodynamic concerns. The design of the propeller blade fillet is a particular 
example illustrated in this research. 

[Nittel 89] presents a method for creating a propeller blade fillet. A number of 
planar sections are specified which are perpendicular to a “nose-to-tail” line. defined 
as a straight line in three-dimensional space connecting the leading edge to the trailing 
edge. The planar sections are also parallel to the y-axis of the coordinate system (see 
Appendix Á for standard propeller geometry definitions). This is somewhat unusual 
since it does not take advantage of the radial nature of the propeller design. The 
complexity of the blade determines the number of sections required. 

Cross-link curves are faired through a small radius fillet and a large radius fillet at 
points of tangency with the hub surface, the blade surface and with each other. The 
cross-link curves form a set of splines which are then merged with the splines repre- 
senting the rest of the propeller blade for use in determining numerically-controlled 
(NC) machine tool paths and cutter instructions. 

[Sabol 83] presents an analytic model for the design of propeller fillets based on 
the standard 3T-T/3 fillet developed at Philadelphia Naval Shipyard. The model 
was developed to provide a complete analytic definition of a blending surface which 
matches as closely as possible the standard 3T-T/3 geometry. A complete analytic 
definition of the fillet which is consistent with the definition of the blade surface 
is required for NC machining. To develop the blade-to-hub fillet, a fillet plane 1s 
established which includes the normal vector to the surface at the intersection of the 
propeller blade with the hub, and a point located along one of the lofting lines of 
the propeller blade (see [Davis 70] for an explanation of the blade lofting procedure). 
The point on the lofting line is called fillet-to-blade intersection. and it is located a 


distance, T, from the intersection of the blade and the hub. The T dimension is the 
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local thickness and is a function of the maximum thickness of the root section and a 
shape factor based on position along the blade-to-hub intersection curve. [he value 
of T is computed for a particular parameter u along the blade-to-hub intersection 
and used to locate the endpoints of the cross-link curve on both the blade and hub 
surfaces. A system of equations is solved to generate a set of points which represent 
the cross-link curves. The equations have four initial conditions which define the 
position and tangent directions at the endpoints of the cross-link curve in the fillet 
plane. The endpoint on the blade is found along the blade lofting line at the distance 
T described above. The endpoint on the hub and the tangent directions are found 
using non-dimensionalized ratios based on the fillet to blade intersection. [he ratios 
were found experimentally to closely approximate the 3T-T/3 fillet. An iterative 
process is used to find the endpoints on the two surfaces and equations are given to 
solve for the tangent directions. Once these four initial conditions are known. the 
system of equations is solved to represent the fillet by point data. The method is 
general in that it makes no requirements on the geometry of the intersecting surfaces. 
It is therefore not restricted to cylindrical or conical hub surfaces. 

The work on feature extraction from B-spline marine propeller representations 
|Patrikalakis 90a] is now mentioned briefly. To implement a 3T-T/3 fillet for this 
thesis. a method to compute the maximum thickness of the hub-to-blade root section 
was required. In [Patrikalakis 90a], a two-dimensional parametric curve is fit to the 
set of points where cylinders coaxial to the propeller axis of rotation intersect the 
blade section. The camber line of the curve is then generated using a system of non- 
linear differential equations. The cylinder selected for use in this thesis corresponds 
to the hub. It is shown that there is a critical point on the camber line and that the 
maximum thickness of the blade section occurs at that critical point. Line segments 
can be constructed from the suction face to the pressure face of the blade. and the 
tangent vectors along the hub-to-blade intersection determined. The critical point 1s 


the point on the camber line where the segment normal to the camber line connects 





points whose tangent vectors are parallel to each other. Taking into account the 
procedure used to construct blade sections and the form ot the thickness function 
which exhibits a single maximum value. the propeller blade section has. in general. 
only one critical point. By finding the critical point on the camber line it is therefore 


possible to compute the maximum thickness of the section. 
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Chapter 3 


Formulation of the Blending 


Problem 


3.1 Introduction 


The blending surface problem is the generation of a smooth transition surface between 
two or more underlving surfaces. The blending surface must match user-defined 
boundary conditions along user-defined boundaries. The user must be able to control 
the interior shape ot the blend. Finally. the surface must be of a form which permits 
interrogation of the surface in and exchange between geometric modelling systems. 


The tasks required to generate a blending surface are: 


l. specification of underlying surfaces 
2. specification of linkage curves 

3. specification of a directional curve 
4. specification of boundary conditions 


. construction of cross-link curves 


Dt 


6. approximation of blending surface 
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i. improvement of approximation 


The specification of the problem will be described in this chapter. The design of the 


blending surface will be demonstrated in Chapter 4. 


3.2 Underlying Surfaces 


As noted in the Introduction (Chapter 1), smooth transitional surfaces are frequently 
required in design problems, whether it be for functional reasons or merely aesthetic 
ones. The utility of an algorithm for the creation of these secondary blending surfaces 
lies. in part, in the types of primary or underlying surfaces the algorithm is able to 
process. Our goal in this work was to specify the underlying surfaces as non-uniform 
rational B-spline (NURBS) surfaces. 

[Hansmann 87] introduced the fundamentals of a general system for the design 
of blending surfaces between primary surfaces of different mathematical representa- 
tion. To demonstrate the method. he used bicubic tensor product B-spline patches 
as examples of free-form surfaces. To reduce computational effort, the B-spline func- 
tions were formulated as third order polynomials instead of using recursive evaluation 
[de Boor 78](see equations 3.6 and 3.7 in Section 3.2.1). A second class of surfaces 
considered were analytically defined quadric surfaces defined in parametric form in 
terms of transcendental functions. For example, an hyperboloid of one sheet with an 
implicit definition of 


F(z.y,2)==+5--1=0 (el) 


Is represented parametrically 1n standard orientation (1.e.. at the origin with principal 


axes aligned with the local coordinate system associated with the surface), by 


ICD COS toD 
lu) = sulu)cosh(v, 
= sv) esI (32) 
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with 0 € u € 2x and =r < v < z. At this point it is worth mentioning the differences 
between implicit and parametric representation. Implicit surfaces represented bv a 


function of the form 


rc 0 (3.3) 


are popular for quadric surfaces. Implicit surfaces have a closed form mathematical 
representation, and as such operations like intersections and surface interrogation can 
be performed exactly. On the other hand, surface point evaluation usually requires 
solution of nonlinear equations. Furthermore. it 1s difficult for the designer to find 
a relationship between the algebraic coefficients which define the surface in implicit 
form and the surface’s geometric properties. [Woodwark87]| notes that high degree 
implicit equations produce surfaces with unpredictable behavior. self-intersections. 
numerical instability, and other related problems. Methods to counteract these prob- 
lems in the representation of free-form curves and surfaces by implicit, piecewise 
algebraic functions, defined in terms of triple product B-splines. can be found in 
[Patrikalakis 88d]. For polynomial (algebraic) surfaces represented in implicit form. 
[Hoffmann 87c| shows implicit polynomial blending surface constructions are possible 
using potential methods. 


Parametric surface representation 1s of the form 


Supe Vr FED Cp esee] (3.4) 


where zx, y, and z are the three Cartesian coordinates of a point on the surface. Points 
on parametric curves and surfaces can be easily computed. Parametric representa- 
tions lend themselves to efficient computations for graphic display and for numerical 
control purposes in NC machining [Faux 79]. The use of vector quantities in para- 
metric equations aids in matrix transformations such as translations and rotations. 
It also allows direct control of the surface. [n the case of blending surface creation. 


parametric surface representations are usually approximations only. Ihe approxi- 


36 





mation 1s controlled by tolerances. Parametric polynomials can represent free-form 
curves and surfaces easily. Since many primary surfaces will be free-form parametric. 
à parametric surface is also a natural choice for the blend. For generality, parametric 
surfaces were selected for the underlying surface and blend representation. 

Non-uniform rational B-spline (NURBS) surface representation was chosen for the 
underlying primary surfaces. The form was chosen because it provides a complete. 
unambiguous mathematical representation able to represent complex sculptured sur- 
faces as well as the classical algebraic surfaces such as quadrics. tori. and surfaces 
of revolution with rational polynomial profiles. Conic sections and quadric surfaces 
can be represented exactly with rational B-splines while retaining the properties of 
integral polynomial B-splines [Tiller 83]. Finally, there exist complete standard data 
formats for NURBS surface representations [Smith 88] which provide for the direct 
exchange of geometrical information between different geometric modelling systems. 
and design and manufacturing systems. 

A brief review of polynomial and rational B-spline theory necessary for this work 


follows. See also [Alourdas 89], [Patrikalakis 89a] and [Tuohy 91b]. 


3.2.1 Polynomial B-Spline Curves and Surfaces 


A spatial, integral B-spline curve has the form 


n—] 
RUDI (3.5) 
=0 


with P the control polygon of n vertices P;, 0 < 2 < n—1, and a knot vector 
MS to ti, te), a set of real numbers such that z; < #;,,,0 <2 < k. The NV; ar(£) 
are scalar valued piecewise polynomials in the parameter t. and are the basis B-spline 


functions evaluated recursively with the Cox-de Boor algorithm (de Boor 72][Cox 72]: 
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l if pe = t = E ; = 
NO = Cay 

O otherwise 
The convention 0/0 = Q is also assumed. Equation (3.5) represents a piecewise. con- 
tinuous polynomial of degree M —1. order M. The curve is continuously differentiable 
up to order .M — 2. If a curve of order .M has no multiple interior knots it is .M — 2 
differentiable. If there are multiple interior knots the continuity order at such knots 


may be reduced. 


The basis functions possess the following properties: 


Noch) > y < l ES n—l (3.51 


nl 
ERE ne en, (3.9) 
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Equations (3.5) and (3.9) form the basıs for the convex hull property. Each interval 
[t;, t;41] of the curve is enclosed within the convex hull of the consecutive .\/ vertices 
affecting the interval. The entire curve is contained within the convex hull of the 


defining polygon vertices. 


The B-spline approximation to a vector-valued, continuous function g(£) in the 


interval [to, £4] is given by 


n— ] 
Vg=) gl&)N:artt) (3.10) 
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The £; are the nodes of the B-spline basis over the knot vector T and are given by 





Si Tp (ett + tiya tiet tima) OSi<n-l (3.11) 


The parametric derivatives of a B-spline curve defined by (3.5) are [de Boor TS) 
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For the same B-spline curve, there exist an infinite number of defining polygons 
with more than the minimum number of vertices. The process of finding an al- 
ternative representation by the addition of new knots to the knot vector is called 
knot refinement or subdivision. Consider a B-spline curve of order M with n control 
polygon vertices. as given by (3.5). The objective is to define the same curve exactly 


with a new knot vector. 


3 d (30, 31, c OD p Sun.) 


obtained from the original knot vector 


Hs = (to, e ets a) 
by the addition of p new knots, (po,....Pp-1) such that m = n + p. The new curve is 
given by 
m-1 
Se) — > opt) (3.14) 
J=0 
where the control polygon vertices C, are such that Ry (t) = Sas(t). A general 


method for knot refinement applicable to any B-spline curve is the Oslo algorithm 


[Cohen 80]. By the algorithm. the new vertices C, are 


n=l Een 
a amc E | (9:15) 
=O O < J L MUS l 


The coefficients ar, are given by the recursion relation 


l a 
NES : (3.16) 
O otherwise 
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e = pecu ba A (3.17) 
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An application of subdivision is the splitting of a B-spline curve at some arbitrary 
point. A new knot can be added at a point with multiplicity M. and the curve can 
be considered to consist of two B-spline curves which. taken together. represent the 
original curve. 

A tensor product B-spline surface is a vector-valued function of two parameters. 
say u.v. The control polyhedron is now topologically a rectangular net of vertices 


eee 2 aml 0a) An l. The two knot vectors 


Da = ee i) (als!) 
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are associated with the parameters u, v. The integral B-spline surface formulation of 


order M in the u direction, .V in the v direction. is given by 


m-in-i 
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By application of equations (3.12) and (3.13) the parametric derivatives of the B-spline 


surface patch with knot vectors (3.18), (3.19) are 
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The primes signify differentiation with respect to u. and the dot with respect to v. 


3.2.2 Rational B-Spline Curves and Surfaces 


Rational B-spline curves and surfaces make use of homogeneous coordinates. For 
example, if a point P = (x,y,z) is defined in 3-dimensional Euclidean space. the 
corresponding homogeneous coordinate is P* — (wz.wy,wz.w), where w is the ho- 
mogeneous coordinate. We assume w > 0. The polynomial formulation of the rational 


B-spline curve in homogeneous space 1s 
n=l 
i—U 


The rational B-spline curve is the projection of the curve defined in (3.26) into 


3-dimensional Euclidean space 


NP B ml) 


R = 
m(t) = Pa een 


Notice that by taking w; = 1 and using (3.9), (3.27) reduces to an integral B-spline 


curve. 
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The generalization to the rational B-spline surface follows: 


au) B .vlv) 
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The homogeneous coordinates w;; are again taken to be positive. As with (3.27). it 
w;; = l the integral B-spline surface patch equation is recovered. 
[Tiller 83] notes that the properties of integral (polynomial) B-spline curves and 


surfaces are easily extended to the rational formulation. 


3.3 Linkage Curves 


The boundary of the blending surface on the underlying surfaces is the linkage curve 
or trim line. The term linkage curve was introduced by [Hansmann 37] and will be 


used in this thesis. There are many ways to specify the linkage curve but most are 


variations of the following: 


e interactively selecting points on the underlying surface (typically by specifying 
them in parameter space and then mapping them to 3-dimensional space), and 
then fitting a curve through these points in the parameter space of the surface 


|Hansmann 37]. This method does not require that the surfaces intersect. 


e specifying points on the underlying surface to simulate rolling ball or spher- 


ical blends and fitting the points with a space curve [Varady 89][Pegna 88] 
[Pegna 90]. 


e specifying an offset curve from the intersection of two surfaces. either along 
geodesic paths on the surface or along isoparameter lines at a constant distance 
[Patrikalakis 89b] [Bardis 39], or as a function of the parameter of the inter- 
section curve [Sabo] 83]. (The geodesic path is locally the path of minimum 


distance between two points on a surface among all paths lying on the surface 





between these points [do Carmo 76]). While the methods are limited by requir- 
ing the surfaces to intersect. the blending of non-intersecting surfaces does not 


often occur in engineering problems. 


Linkage curves which are not isoparameter lines will generally be of high degree. 
Depending on the method used to define the curves. the linkage curves will be either 
high degree B-splines, procedural curves, or curves defined in the parameter space 
of the underlying surfaces. No matter what method is used to specify the linkage 
curve, care must be taken to ensure the curve does. in fact, lie on the surface with 
sufficiently high accuracy. 

The linkage curves must be properly parameterized. When constructing the cross- 
link curves it is necessary to establish some type of correspondence between the linkage 
curves. The usual method for establishing this is by a correspondence of equal param- 
eter values. If the linkage curves are not carefully parameterized, however, there is 
a possibility that this algebraic isoparametric correspondence will link inappropriate 
points causing undesirable twist and intersections [Koparkar 91] [Pegna 89]. Inter- 
estingly, [Koparkar 91] suggests that the parameterization of the underlying surfaces 
and the linkage curves is unimportant. A suitable blend requires the proper spatial 
relationship between the curves and the underlying surfaces. To ensure proper corre- 
spondence, he uses Euclidean projections from the intersection curve of the normal or 


fanout surfaces onto the underlying surfaces to establish a geometric correspondence 


between the linkage curves. 


3.4 Directional Curves 


Most of the blending surface design methods construct the blend surface with a family 
of cross-link curves. The cross-link curves usually connect the linkage curves through 
a specified correspondence scheme and are typically designed as planar curves. The 


use of a directional curve is one method of defining the orientation of this planar 
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curve. The directional curve may be defined as: 


e the intersection curve of the surfaces to be joined, 
e the intersection of offset surfaces or fanout surfaces of the underlying surfaces. 


e a free-form space curve specified by the designer to achieve the desired effect 


both at the boundary and interior to the blending surface. 


e the orthogonal projection of one linkage curve onto the other surface. 


The usual means of using the directional curve is through an algebraic corre- 
spondence of parameters between the linkage curves and the directional curve. The 
technique essentially determines the direction of tangent vectors to the surface at the 
linkage point. Alternative techniques are to use orthogonal projections of the direc- 
tional curve (in this case the directrix) onto the underlying surface, or to establish 
a geometric correspondence between a Euclidean projection of the directional curve 
onto the surface. 

An alternative to a directional curve is to define a function which determines the 
direction of the tangent vectors at the linkage point. [Hansmann 87] requires the 
cross-link curves to be orthogonal to the linkage curve. In an automatic marine pro- 
peller fillet example demonstrated in Chapter 5. the cross-link curves are constructed 
in “fillet planes.” Each fillet plane contains a vector normal to the propeller blade 
surface at the intersection of the blade and hub. and a point on an isoparameter 


curve extending along the blade surface from the hub to the tip (see Section B.2). No 


directional curve is required. 


3.5 Boundary Conditions 


Boundary conditions specify the degree of continuity between the blend surface and 
the underlying surface along the linkage curve. The simplest type of boundary con- 


dition is position continuity. This simply ensures the edge of the blend surface does 
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in fact lie on the underlving surface along the linkage curve. to within the specified 


tolerance. No further conditions are required of the blend. 


3.9.1 Tangent Plane Continuity 


Most blending algorithms require. beyond positional. tangent plane continuity across 
the linkage curve. [Pegna 89a] refers to this as the first order blending problem. To 
satisfy tangent plane continuity, the tangent planes to the blend and the underlying 
surface must match at all points along the linkage curve. This condition is satisfied 
during construction of the cross-link curves. It may be automatic such as in the case 
of circular or spherical sweeps. where the generatrix contacts the surface at a point 
of tangency. More general methods determine the tangent plane to the underlying 
surface at a point on the linkage curve, and construct the cross-link curve by initiating 
the linkage curve in the tangent plane in the direction specified by the directional 


curve. A bias factor can be applied to control the effect of the tangency condition 


into the interior of the blend. 


3.5.2 Curvature Continuity 


Curvature continuity, or second order blending, is achieved when normal curvatures 
agree in every direction at every point along the linkage curve. Before proceeding 


with the requirements for curvature continuity, a review of some differential geometry 


concepts is presented. 


3.5.2.1 Review of Differential Geometry and Curvature 


We take the curve u = u(t) on a parametric surface r 2 r(u.v) where 


r=r(u.v)=| y(u.v) (3.29) 





represents one of the underlying surfaces in the blending problem. and 
BERN = (3.30) 


the linkage curve on that surface. A specific point on the linkage curve in three 


dimensions will be referred to as r(t). 


The unit vector normal to the surface is given by 


Ey A Ty) 


= c 
IU rr. | Mors 


where ^ denotes the vector product. A unit tangent vector to the linkage curve. r’. 


is given by [Faux 79] 


n UE A (308, 
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Recalling that the arc length of a parametric curve is 


t 
s= e (3.34) 
to 
the squared arc element is 
LX O o dr (3.39) 
or. more simply. 
A A at (3.36) 
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where 
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are the elements of the first fundamental form. This is also referred to as the first 
fundamental matrix or the metric tensor. It ıs a measure of the distance in three 
dimensional space that exists between neighboring points on a surface. 


The curvature of a curve parameterized bv arc length is 


Exe (3.40) 
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The radius of curvature is 
l 
K 


Using (3.32) for the curve u(t) = [u(t), v(t)| on the surface r = r{u.v), we obtain 


Hn 2 R S oD e E : 
D reduc PLU 8 (3.43) 


Taking the component of (3.43) in the direction of the surface normal n. 
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and recalling that n-r, =n-r, — 0. we find the second fundamental form or tensor 


Dosen c UMEN E ND TIS 
where 
AO EA Pus (3.46) 
iM — A | a Tr, (Sala 
O SA Pos (3.48) 


The second fundamental tensor provides a measure of the curvature of a surface at 
a point. It can be shown that the normal curvature &4 of the surface at a point r(t) 
along the linkage curve in a given tangent direction u,v is simply the ratio of the 


second fundamental form to the first fundamental form. 


4 u +2Müo+ NÜ 


Z- Te 3.49 
“Ew +2Fu0 + Ge ae 


A normal section is the plane containing both a unit vector t = [ù, ù] emanating from 
a point r(t) in the tangent plane. and the normal vector to the surface at that point. 
Equation (3.49) therefore gives the normal curvature for any normal section at a 
point r(t). [Faux 79] presents the notion of signed curvature for planar curves and for 
surfaces. For planar curves. the sign of curvature is chosen as required for the applica- 
tion. For surfaces. the curvature is assumed to be positive if the curve of intersection 
of the normal section with the surface has its center of curvature in the direction of 
the positive surface normal vector. The principal curvatures are the maximum and 
minimum values of (3.19). and are denoted x1,k2. The product X = x&4&» is the 
Gaussian curvature, and the mean curvature is H — s(&1 +k). The directions cor- 
responding to the principal curvatures are called the principal curvature directions. 


A line of principal curvature of the surface is a curve on the surface such that. at 
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all points r(t) along the curve. the tangent line of the curve 1s a principal direc- 
tion. It can be shown that the lines of principal curvature are orthogonal [Farin 83] 
ido Carmo 76]. 

An asvmptotic direction of the surface r(u,v) at a point r(t) is a direction in the 
tangent plane for which the normal curvature is zero. A geometric interpretation of 
asymptotic direction is the equation of a conic section known as Dupin's Indicatrix. 
Knowing the principal curvatures. the normal curvature A, of the section formed by 


the unit tangent vector t and the surface normal vector is 
Ka =k, cos B -- &.sin^ Ü (3209 


where 0 1s the angle between the unit tangent vector t and the direction of the 
maximum principal curvature. Equation (3.50) is the classical Euler's Formula, ex- 
pressed in terms of the orthonormal basis formed bv the principal curvature direc- 
tions. Dupin's Indicatrix is the union of conic sections in the tangent plane such that 
Ky = +1. When &, and &5 have the same sign at r(t), the Gaussian curvature is pos- 
itive, r(t) is an elliptic point, and Dupin’s Indicatrix is the graph of an ellipse. If &i 
and &2 have opposite signs, the point 1s hyperbolic. and Dupin`s Indicatrix gives two 
hvperbolas with a common set of asymptotic lines. If one of the principal curvatures 
is zero. the asymptotic lines are parallel straight lines and the point is parabolic. 


For an extensive treatment of these concepts refer to [do Carmo 76] and [Farin 88]. 


3.5.2.2 Conditions for Curvature Continuity 


With these concepts in mind. we now consider the problem of curvature continuity of 
adjoining surfaces. Pegna and Wolter [Pegna 89a] present the following formulations 
for curvature continultv: 

Ihe osculating paraboloids of the underlving surface and the blending surface 
must match at all points along the linkage curve. The osculating paraboloid for 


a surface S with the ry plane as the tangent plane at a point p on S is given bv 


+9 


o 





do Carmo 76]: 


l : 
z= (a fee + 22Y fe, ty’ fy) (cron) 


where = = f(z. y) 1s the equation of the surface. 


If the blend juncture is already shown to be tangent plane continuous, the following 


criteria may be used: 


l. The Dupin's indicatrices and curvature signs of both surfaces must agree along 


the linkage curve. 


2. The asymptotic directions of the two surfaces agree at all points along the 


linkage curve. 


J. The second fundamental tensors. formed as basis independent. multilinear maps. 


agree at all points along the linkage curve. 


i. The principal curvature directions and principal radii of curvature for the sur- 


faces agree along the linkage curve. 


Finally, as noted in Chapter 2, [Pegna 89a] shows that two surfaces which are 
tangent along a first order smooth linkage curve are curvature continuous if their nor- 
mal curvature agrees in a direction other than the tangent to the linkage curve for all 
points along the linkage curve. This formulation is called the Linkage Curve Theorem. 


Since it is important to our work, this theorem is stated formally from [Pegna 89al: 


Theorem 3.1 Given two C*-continuous surfaces S4 and S», and a Cl-continuous 
linkage curve E(t) such that E' Z 0 for allt. t € T with T C R. where X is the set of 
real numbers. The linkage curve represents the tangent plane continuous intersection 
of the two surfaces. If at some point E(tg) € E(T) for some tangent direction X 
linearly independent from the tangent direction E'(tg) the normal curvatures of S, 


and S, agree, then the normal curvatures of Sı and S, agree for all tangent directions 


at the point E(tg). 





Chapter 4 


Blending Surface Generation 


4.1 Introduction 


The method for blending surface generation implemented in this thesis approximates 
a smooth transitional surface as a B-spline surface with user-specified boundary con- 


ditions across the linkage curves. The general method is as follows: 


l. Approximate the linkage curves with integral cubic B-spline curves. The knot 
vectors of the approximated linkage curves are then merged to form a basis for 


the design of the cross-link curves. 


Lo 


Construct the family of cross-link curves based on the type of blending surface 
desired and the specified boundary conditions. linkage and directional curves. 
and bias functions. The parameter values of the cross-link curves are the nodes 


of the merged knot vector of the linkage curves. 


3. Create the approximating blending surface by lofting the cross-link curves over 
the merged knot vector of the approximating linkage curves. [Tuohy 91b] shows 
that by using the merged knot vector and the nodes of the linkage curves. the 


boundary curves of the blend will be identical to the curves approximating the 


linkage curves. 


1. Evaluate the approximation by checking continuity errors with adaptive sam- 
pling of the blending surface and underlying surfaces. The degree of sampling 


is dependent on the blend continuity specified across the linkage curves. 


5. Improve the approximation by the addition of cross-link curves based on the 


adaptive sampling results until the error 1s within the specified tolerances. 


4.2 Linkage Curve Approximation 


The linkage curve defines the extent of the blend on the underlying surfaces. It 
can be expressed as a B-spline curve, a curve defined in the parameter space of the 
underlying surface. or as a procedural curve. No matter how the linkage curve is 
specified, it must in fact lie on the surface to ensure the blending surface joins the 
underlying surfaces with proper position continuity. We use a method presented bv 
[Wolter 91] and [Tuohy 91b] to approximate the linkage curve with a cubic integral B- 
spline curve to within position. tangent angle and relative curvature tolerances when 
compared with the original linkage curve. We assume the linkage curve 1s expressed 
as a vector-valued function, g(t), and seek an approximating integral B-spline curve 


as in equation (3.5), 


n—] 
Ra(t) 2 V P,B,u(1) EN 
=U 


which interpolates the function g(t) at the nodes £; of the knot vector 


T = Bot ug), (Gee 


normalized such that the overall parameter interval spans [0.1]. Selecting the nodes 


Eun (3.11). 


l 
G= Woy itt LL t. (4.3) 





the vertices P; of the control polvgon of the approximating B-spline may be computed 


by solving the system 
n= 1 
ee Ee SS n— 1. (4.4) 


J=0 


which can be written in matrix form as 


[NIP] = {g} (4.5) 
where 
[NV] E p oe a (4.6) 
ea ce) Cum 
(P) 2 [P5....P. i]! (4.8) 
elle ecl (4.9) 


[de Boor 78] shows that the matrix [| V] is non-singular and therefore invertible since 
we evaluate the system at the nodes of a monotonically increasing knot vector. The 
system of (4.5) can thus be solved for the control polygon using standard linear algebra 
techniques. 

The initial approximation is then improved by adaptive sampling. The method is 
explained here because an adaptation of the method is used to improve the blending 
surface. The approximation is tested at midspan between the nodes for position error 
using local error bounds. The approximation method ensures the position error at 
the nodes is zero. [Wolter 91] and [Tuohy 91b] show that the maximum position error 
occurs near the midspan between two adjacent nodes. The error function. therefore. 
is simply 


R Ese mes E NEM dene >. CEU 


cm 





[f the position test is satisfied. the angle o is formed by the tangent vector to the 


approximating B-spline and the tangent to the linkage curve. 


arccos E j ul) 
JRA, (us) ([g/(u:)| 


En 
< €2, u = E for 


A 


u 








where 7; is defined as in (4.10). In practice. it is usually convenient to test the angle 
difference at the nodes as u; = €; as the distance between successive knots decreases 
during adaptive sampling. 

In similar fashion. relative curvature at the nodes of the approximating curve 


Ir’ A r| 


K = rp 


is tested such that 
Kgl(€;) — KR(&) 
Kg (€i) 


For both first and second order tests, the use of geometric properties such as angle 


m (4.12) 








or curvature of the approximation and original curve is preferred [Tuohy 91b] since 
such properties are independent of the curve parameterization. At points where the 
error test fails, t, is added to the knot vector 7. With this enriched knot vector. a 
new approximating B-spline curve is constructed. The procedure is repeated until 
the approximating curve passes all error tests. 

From the above discussion, it is noted that improvement of the approximation 
is restricted solely to knot vector refinement. The shape of any B-spline curve is 
derived from its control points and the basis functions. The knot vector has a sig- 
nificant influence on the basis functions [Rogers 76] for a given curve order. Since 
the approximation scheme interpolates the linkage curve at the nodes exactly, knot 
vector refinement offers the best means of improving the approximation without loos- 
ing position accuracy. Furthermore. prudent selection of the initial knot vector will 
enhance the efficiency of the approximation. Several methods for selection of initial 


knot vector are available [Hartley 80] [Hoelzle 33] [Hoschek 87]. 
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4.3 Construction of Cross-Link Curves 


Two techniques were examined for the construction of cross-link curves. The first 
algorithm 1s a general one which creates a family of cross-link curves derived from 
the linkage curves. the boundary conditions across the linkage curves. a directional 
curve. and a bias function. The second algorithm is actually a special case of the 
first, where the linkage curves and cross-link curves are constructed automatically 
to satisfy specified geometric conditions. The second algorithm 1s presented as an 
example in Chapter 5. 


The cross-link curves are integral B-spline curves of the form 


C(t) — >> Ba), er 


il 


where M is the order and m the number of vertices of the control polygon. P. The 
cross-link curves are constructed at the nodes of the common knot vector created by 
merging the knot vectors of the approximations to the linkage curves. As noted at the 
beginning of this Chapter, this procedure ensures the boundary curves of the lofted 
blending surface are identical to the approximating linkage curves. The boundary 
conditions are used to determine the order of the cross-link curve. the number of 
control points, and the knot vector. Nine combinations of position. tangent plane 
and curvature continuity are possible. 

The simplest case is position continuity on both underlying surfaces. The cross- 
link curves are of second order and have two control points with a knot vector of 
(0.0.1.1). The control points are simply the points on the linkage curves associated 
with the particular node of the common knot vector being evaluated. This type of 
cross-link curve is useful as a preliminary check of the linkage curve correspondence. 

The next two cases require position continuity on one surface. and position and 


tangent plane continuity on the other surface. A third order B-spline curve with three 


Wt 
Qt 





control points. 


C(t) = 37 PiNalt), (4.14) 
¡=0 


is used with the knot vector (0.0.0.1. 1,1). These curves are all of Bézier form. a 
special case of B-spline, and can be raised to cubic B-splines easily |Cohen 35]. 

The remaining combinations use integral cubic B-splines (.M = +). One case 
requires position and tangent plane continuity on both surfaces. Two other cases 
require position continuity on one surface, and position. tangent plane and curvature 


continuity on the other. The cross-link curve for these cases is of the form 
3 
Ce o (4.15) 
p 


with knot vector (0.0,0.0. 1.1.1, 1), and is again of Bézier form. The three cases 
differ from one another only in the manner in which the control points are computed. 
Two cases require up to tangent plane continuity on one surface, and up to cur- 


vature continuity on the other surface. For these cases a fifth control point is added. 
= : 
C= DI ade (4.16) 
= 


with knot vector (0,0.0.0. 5, 71 120200500 i, 1,1.1,1). For convenience the 
cross-link curve C(t) is parameterized such that C(0) interpolates surface S, with 
linkage curve E,(t), and C(1) interpolates the second surface S4, with linkage curve 
F(t). The knot vector is chosen such that the internal knot is closest to the surface 


with the curvature continuity boundary condition. 


The final case requires up to curvature continuity conditions with both underlving 


surfaces. The cross-link curves now take the form 
Gt — s pU. (UNES 
m 


with six control points and knot vector (0,0, 0.0. = 5, 1.1.1.1). Having selected the 
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proper form of the cross-link curve required to satisfy the boundary conditions. the 


Mirves are constructed. 


4.3.1 Position Continuity 


Position continuity is guaranteed to a specified tolerance by evaluating the linkage 
curves at the nodal points of the common knot vector. We label these points on the 
linkage curves E;(£;), J — 1,2, corresponding to surfaces S;. A common knot vector 
T. derived from the approximation of the linkage curves, is evaluated at the nodes £i, 


as defined in (3.11). The cross-link curve associated with the node £; has the form 


C(t) = e P;N, (t) (4.18) 


2=0 


To satisfy position continuity, then, we require that 


DN SE E 


Pn-ı - E) (4.20) 


4.3.2  Tangent Plane Continuity 


langent plane continuity for each cross-link curve requires two control points asso- 
ciated with each boundary. The first establishes position continuity. The second 
control point is located ın the plane tangent to the primary surface so that the cross- 
link curve 1s tangent to the surface at the linkage curve. The position of the second 
control point in the tangent plane is determined by the directional curve and a scalar 
bias function. This function is used to control the magnitude of the tangent vector 
of the cross-link curve and thus influence its shape. 

The direction of the vector in the tangent plane of the primary surface at E;(£;) 
is based on a technique developed by [Tuohy 91b]. The method is illustrated in 


Figure 4.1. The direction t; lies along the intersection of the plane tangent to the 


= 
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Figure 4.1: Tangent plane construction 


primary surface at E,(€;), j= 1,2, with a plane containing corresponding points on 
both linkage curves and the directional curve. The direction is found by computing 
the cross product of the normal vectors to the two planes. 

We first compute the unit normal vector to the plane formed by the points EE p 
where E, and E; are corresponding points on the linkage curves and D is the assocı- 


ated point on the directional curve. The unit normal vector 1s computed as 


(E; —- Ei) ^ (D — E) ( 
(E: — E1) A (D — E; )| 


no = 


4.21) 


The unit normal vector n; to the underlying surface S;, ; — 1,2 1s computed as 


SIN WS. 


n = —————. 
IS. AS,| 


In the event one of the parametric derivatives of the surface is zero, (1.e.. a degenerate 
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patch). the unit normal vector to the surface is found by either 


Sey 
es an fS,=0. 4.23 
n; IS., ^ S, i or ( ) 
Su A Sus » 
> IS. T2 
"i Beasa an 


as suggested by [Faux 79]. 


The control points of the cross-link curve required to establish tangent plane 


continuity at its end point on the linkage curve E,(t) are, therefore. 


Pee Bie) (4.25) 
DUE Poz ti (4.26) 
where tı = q(mAnı) (4.27) 


where qi is a positive bias factor which can depend on £;. The control points for the 


other end of the cross-link curve associated, with Ez, are 


eee = CES (4.28) 
lc TP (EU (129) 
where t = q(no^ n3) (4.30) 


where q2 1s also a positive bias factor which can depend on £;, and be defined inde- 
pendently from qı. The + in equations (4.26) and (4.29) provide a way to ensure the 
tangent vectors are oriented in the proper direction. Depending upon the order in 
which the cross products are taken and the direction of the surface normal vector, the 
tangent vector will point either towards the directional curve or away from 1t while 
still lving along the intersection of the tangent and surface normal planes (Figure 4.2). 
Without a directional check. cross-link curves and blending surfaces may arise such as 


shown in Figure 4.3. To account for this we compute the angle between the tangent 


ST 
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Figure 4.2: Possible directions of t 
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Figure 4.3: Blend with wrong tangent direction 


vector t; and the vector D — E;, 


op) | 
Qj — arccos ————— ——-—. COM 
t, |D n E;| 
and the angle between the opposite tangent vector. —t;, and D — E;. 
—t;-(D—E,) 
O = arccos —————, (4.32) 
t;||D — Ej, 


If o, < 01, the tangent vector —t, is selected, thus determining the sign in equa- 
tions (4.26) and (4.29) and ensuring the correct tangent vector is chosen. The benefit 
of this technique is that the choice of proper tangent direction is automatic and the 
tangent is guaranteed to point towards the directional curve. Thus. if a blending 
surface such as that shown in Figure 4.3 is desired. all that is required is that the 
directional curve be located to select that tangent orientation. 


The factor q; in equations (4.27) and (4.30) represent the scalar bias functions. 
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The bias function controls the shape of the cross-link curve by controlling the mag- 
nitude of the tangent vector. It can be a constant value or a function q;(t) which 
varies with the parameterization of the linkage curves or directional curve. Examples 
of constant and variable bias functions are given in Section 5.2. 

The method presented above (equations (4.25) to (4.30)) allow the definition of 
planar cross-link curves. each of which is in a plane E, ED controlled by the paramet- 
ric correspondence of the linkage and directional curves. This method is at variance 
with the method developed by [Hansmann 37], who chooses vectors t; and t» to be 
orthogonal to the linkage curves (as well as to n; and n5) and does not. therefore. 
have planar cross-link curves even for up to tangent plane continuous blends. Our 
experimentation indicates that the construction of planar cross-link curves with tan- 
gent vectors at the ends not exclusively controlled bv the linkage curves but rather 
from the directional curve provides easy control of the shape of the blend and good 


flexibility in the design of complex free-form blends. 


4.3.3 Curvature Continuity 


To satisfy curvature continuity boundary conditions, we rely on differential geometry 

properties of the underlying surfaces along the linkage curves (see Section 3.5.2.1). 
The cross-link curve is of the form given by equation (4.13). Evaluating this 

equation and the first and second parametric derivatives given in equations (3.12) 


and (3.13), we get the following system of equations along E(t): 


Po = Ede) (2.99) 
i ¢ 
P, = Po+C; (0) (==) (4.34) 
len AA => b — tł ru to 
P, = P, n (o) ( + M-2 LM 1 | 
(tar — to) (tater — te) ae 
O uis d 1.3 
gi(0) ( (M — 1M —2) a 





For control points associated with the second surface and Ex(t). 


Pm-1 = Ex) (4.36) 
Im M-2 ae = 

De S= Paas QU ——— (4.37) 
Im M-— Im M—3 lm- p Im-: 

O ene fim fe). 


^ (1) d (4.38) 


i (M —1)(M —2) 
Pesbeiore, E.(&,), ) = 1,2 represent the linkage curves evaluated at the nodes £; of 
the common knot vector. 5, of the approximations to the linkage curves. and C is 
the corresponding cross-link curve and its parametric derivatives. The subscript É; 
denotes association with a particular node. The knot vector for the cross-link curve 1s 
T — (tg,t4,... t LM 1) 88 specified in Section 4.2 above. The control points Po, P, 
and P4.,,P,,.? are found using the method shown in Section 1.3.2 to satisfy position 
and tangent plane continuity. From equations (4.33), (4.34), (4.36) and (4.37), we 
can therefore solve for C; (0) and C; (1): 
LO) = (Py = Po) — (4.39) 
E ES Det’. (4.40) 


Du m) = Dd | | 


To find the remaining control points P; and P4. 5, we need information concern- 
ing second derivative of the cross-link curve at its end point. We find this using a 
generalization of the method presented by [Hansmann 87]. Specifically, we obtain the 
second derivative of the cross-link curve, evaluated at the end points E,(£;), from the 
curvature of the normal section of the underlying surface which contains the cross-link 
curve tangent vector and the surface normal vector at E;(£;). 


We first find the curvature of the normal section curve using equation (3.49) 


— Li? Mui Ni? 


oN) = Fo cr 


(4.41) 
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where E. F.G. L..M,.N are the coefficients of the first and second fundamental forms 
defined in equations (3.37) through (3.39), and (3.46) through (3.48), and (à, v) the 
direction of the normal section curve in the parameter space of the primary surface. 
The dot signifies differentiation with respect to the parameter t of the curve on the 
surface formed by the intersection of the normal section plane and the underlving 


surface. The underlying surface, of the form 


m—1n-l 


Su N(u,v) = $^ 5, PijBiu(u)B; (v) (4.42) 
i=0 7=0 
is evaluated for position. first, second and cross-derivatives recursively using the Cox- 
de Boor algorithm and equations (3.21) through (3.25). The results are used to solve 
for the first and second fundamental coefficients. 

To find the required tangent direction (u,v) of the normal section in terms of 
the surface parameters we recall from differential geometry (Figure +.4 and equa- 
tion (3.32)) that a nonzero vector t is tangent to the surface S at the point E(€,) if 
there is a regular curve r — r(t) on the surface through E(£;) such that t — r. If 


r — r(t) 2 r(u(t), v(t)), then the tangent vector to the surface is 


t=r,-utr,-v (4.43) 


The vectors r, and r, are the first parametric derivatives of the surface at E(£;) 
and are known. We also know t from the solution to the tangent plane continuity 
problem in Section 4.3.2. Denoting the components of the vectors as t = (tx,ty,tz). 


IE Ya -u) andr,=(r,y.,-,) in the global coordinate system. equation (4.13) 


can be written in matrix form as 


>< 

E 
e 

E 
= 


ty E u Uv (4.44) 
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Figure 4.4: Tangent vector in terms of surface parameters 


This is a linear system of the form Ar = b. We solve for [ù ù]! using a least squares 
technique [NAG 89] by taking 


t= e, (0) (4.45) 
and 


t = C4 (1). (4.46) 


We then evaluate (4.41) to obtain the normal curvature &4(€;) of the primary surfaces 
along the normal section formed by fù, ù]. 

Having obtained the normal curvature of the primary surfaces at the points E, (£;) 
and directions t;,, (3 = 1,2), we now proceed to construct cross-link curves which 
have the same curvature at these points and which are also locally normal sections of 
the blending surface to be generated. Because the tangent plane continuity 1s already 
enforced, the relevant normal section planes of the blend at E, and E; are defined by 
(t1, n1) and (t5, n3), respectively. 

]o relate the normal curvature to second derivative of the cross-link curves at their 
end points, we consider a local orthogonal coordinate system with the origin at E,(£;). 
the z-axis parallel to the direction of the tangent vector t;, the y-axis parallel to the 


primary surface normal vector n;, and the z-axis defined by t; ^n;. We now compute 
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the components of the cross-link curve second derivatives at their end points in these 
two local coordinate systems. We use equation (3.11) for r 2 C&4(t) 2 izi(t). yi(t).0 


for t near 0 or 1. In this wav. we obtain the relation 
za — v2; | = ®alCz,(r)|° (d) 


for rT — 0 or 1. Using equations (4.45) and (4.46) and decomposing t along the local 
axes so that t = [tz,t,,t.] where t. = 0 and |t| = (¢2 + ¢?)!/?. then equation (4.47) 


becomes 


yr = | = ER rS 


Equation (4.43) is satisfied bv choosing 


fi 
Tı 


tl ts (4.50) 


ff 
Y 


These components of second derivative [z/. yj, 0) are then transformed to C7 (0) or 
C; (1), as appropriate, in the global coordinate system using the direction cosines of 
the axes of the local coordinate system. The final control points P2 and P,,-3 are 
then solved for using equations (4.35) and (4.38). In this manner. because Cz (0) 
and C,,(1) are vectors in the planes (t,,n, and (t2,n2) the resulting control points 
will no longer necessarily be contained in the plane E,E;D. Therefore. for curvature 
continuous blends the resulting cross-link curves are non-planar. 

In the implementation of this work described in Chapter 5 and in order to main- 
tain planarity of cross-link curves and their easv control. we developed an alternate 
method. In this method we impose the curvature values at the ends of the cross-link 
curves as obtained from normal sections of the primary surfaces in the directions t, 
and t2. However. we use local coordinate systems containing the plane E ED Mes 


local z-axis in the t; direction and local y-axis in the no A t; direction). In this man- 
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ner. the resulting cross-link curve is planar in its entirety and contained 1n the plane 
E, E-D. However. this curve is no longer locally a normal section of the blend at E, 
or E, unless no: n, = 0 or no. n; — 0 and. therefore. not all conditions of surface cur- 
vature continuity are generally satisfied. This implementation of Chapter 5 resulted 
in easily controllable blends of smooth form which can be characterized as providing 


continuitv of curvature of cross-link curves. 


4.4 Approximation of the Blending Surface 


The blending surface is created using a lofting technique where the cross-link curves 
are lofted using the common knot vector from the merging of the knot vectors ot 
the approximations to the linkage curves. The cross-link curves become 1soparameter 
curves of the blending surface. As noted earlier. use of the common knot vector and 
cross-link curves at the nodes of the knot vector ensures the boundary curves of the 
blending surface will be identical with the approximations of the linkage curves. 

A technique similar to the improvement of the approximation to the linkage curves 
1s used to improve the accuracy of the blending surface. The surface is sampled at 
the midspan between the nodes of the common knot vector and checked against 
tolerances to the specified boundary conditions [Tiller 83]. Position. tangent and 
curvature continuity along the linkage curve is guaranteed by the lofting technique. 
To check for tangent plane continuity of the blend, we compute the angle between 
the unit normal vector to the blending surface, n,,, along the boundary of the blend 


with the unit normal vector to the associated underlving surface. M 


HH. : 
O = arccos (ee) (4.51) 


In, In, | 


with 7 = 1,2. to indicate the underlying surface and linkage curve and b; the appro- 
priate boundary curve of the blending surface. Sampling is at the midspan between 


the nodes as defined in (4.10). If the angle found in (4.51) is greater than the specified 
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tolerance. «1, the value t; is added to the common knot vector. 

To test for curvature continuity, we extract an isoparameter curve in the cross- 
link direction from the blending surface. again at the midspan of the nodal points. 
The isoparameter curve comes from evaluating the surface over the complete range of 
one parameter maintaining the other parameter constant. The blending surface was 
constructed such that the cross-link curves become isoparameter curves of the lofted 


blend surface. Let the blending surface be given by, 


m-l] n=l 


Bee) — c UPS B(u)B;v(v). (a) 


i0 321) 


where u.v are the parametric directions along the linkage and cross-link curves. The 


Isoparameter curve extracted from 1t is of the form 


n=l 
I(v) = m Bao (1.53) 
J=0 
with 
M-1 
R; = > O Bu) [m 
1=0 


We then compute the curvature of the isoparameter curve at the endpoint with a 


curvature continuitv boundarv condition. 


SEE ATE) 


= 155 
AI Top (4.55) 


For comparison. we construct a cross-link curve at the same midspan nodal value and 


compute its curvature, x.ı, in the same manner. The curvatures are then compared 


to a relative curvature tolerance, e,, defined as 


I&r p Kell 


ler] + [Keil 


As before, if the test fails the midspan point is added to the knot vector. 
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The improvement process for curvature continuity compares the blending surface 
at a point with the cross-link curve which would be added to that point 1f required. 
If the new cross-link curve does not improve the surface. it 1s not added. The im- 
provement of the blend surface tests tangent plane continuity first if tangent plane or 
curvature continuity 1s specified. If the tangent plane test fails the point is added to 
the knot vector and the next interval is sampled. Curvature continuity is only tested 
if tangent plane continuity is satisfied first. Once all intervals are tested the blending 
surface 1s reconstructed with the new knot vector. The process is continued until all 
continuity tests are satisfied. 

As a final remark. we assert that a blending surface lofted from cross-link curves 
constructed as presented in Section 4.3 while also being normal sections of the blend 
at E, and Ez, is curvature continuous by the Linkage Curve Theorem (Theorem 3.1). 
The tangent directions at the end points of the cross-link curve are guaranteed to be 


independent from the tangent direction of the linkage curve for appropriately oriented 


linkage and directional curves. 


69 





Chapter 5 


Numerical Examples 


We now illustrate the generation of blending surfaces with several numerical ex- 
amples. The first examples demonstrate just position continuity along the linkage 
curves. Some examples are given to illustrate some important concepts concerning 
blending surface shape control with the bias functions. These examples use simple pla- 
nar and quadratic integral and rational surfaces which permit a clear demonstration 
of the results. The last sequence of examples illustrate tangent plane and curvature 
continuous blending surfaces with complex surfaces representing marine propellers. 

Ihe routines to compute. interrogate. and displav the blending surfaces were 
developed for use in the MIT Ocean Engineering Design Laboratory geometric mod- 
eling system, Praziteles. This system provides an environment for both graphical and 
computational interaction between a designer and the Silicon Graphics IRIS 4D 120 
GTXB workstation. The basic structure of the system is described in [Hottel 89]. 
with more advanced versions documented in [Alourdas 90], [Patrikalakis 90a] and 
[Hottel 91]. The routines are written in the C programming language with calls to 


FORTRAN subroutines in the NAG Mark 13 FORTRAN library [NAG 89] where 


required. 


5.1 Position Continuous Blends 


Position continuity is the simplest type of blend. The term “blend” is used loosely 
since the surface generated does not provide a gradual transition between the un- 
derlying surfaces but simply connects them. The primary advantage of the position 
continuous blending surface is that it is created very quickly. In our implementation. 
a preliminary position continuous blending surface is constructed after specification 
and approximation of the linkage curves to verifv proper parameterization of the link- 
age curves. If the linkage curves are designed independently there is à possibility of 
introducing unwanted twist or self-intersections into the blending surface. This is 
illustrated in Figure 5.1. The first surface is an integral B-spline planar patch of the 


form given in equation (3.20), with knot vectors in the u and v directions of 
2 (000.1. 1) VE OU 
and control polygon vertices in three-dimensional space 
IS OO EP) =(0,1,0) Pii = (1,1,0) 


The linkage curve is linear with knot vector and control polygon vertices in the pa- 


rameter space of the patch 
Een 0.2 0.00 P, = (0.7.1.0) 


Ihe second surface patch ıs also planar with the same u.v knot vectors and control 


polygon vertices of 


Poo — OO) Po; zm (0, 0. 0) Pio = Lu. 1) Pii = (0.0. Iu 





à DA À 
SO 


Figure 5.1: Poorly parameterized linkage curves 


The linkage curve is linear with knot vector and control polygon vertices in the pa- 


rameter space of the patch 
Eu UNUS Do — Oro. 0.0) FP, = (0.5, 1.0) 


The directional curve is the intersection of the two surfaces expressed in three- 


dimensional space as a linear spline with knot vector and control polygon vertices 
ES oS = 0.0.0) D, = (1.0.0) 


À minimum number of cross-link curves with position continuity boundary conditions 
along both linkage curves are constructed. Verification of appropriate correspondence 
of the linkage curve parameterization is a visual check. As shown in Figure 5.1. 
the linkage curves are parameterized in opposite directions. introducing a twist in 


the blending surface. To remove the twist. the linkage curve on the second patch 





Figure 5.2: Satisfactory parameterization check 


is reparameterized with new control points Po = (0.5.1.0), Pı = (0.5.0.0). and the 


visual correspondence check performed again. The results of the new parameterization 


are shown in Figure 3.2. 


5.2 Bias Function Effects in Tangent and 
Cross-Link Curve Curvature 


Continuous Blends 


As noted in Section 4.3.2. the bias function is used to control the shape of the blending 
surface, and can be given as a constant value or as a function based on the param- 
eterization of the linkage curves or directional curve. The bias function determines 
the magnitude of the first derivatives of the cross-link curves at their endpoints. À 
constant value bias function was implemented by [Bardis $9]. The technique requires 


the distance between corresponding points on the linkage curves to be relatively con- 


13 








Figure 5.3: Constant scalar bias function, q; — 0.4 


stant. As shown in Figure 5.3. if the distance varies significantlv, a constant bias 
factor. suitable where the linkage curves are far apart, can result in self-intersections 
or cusps in the cross-link curves and the resulting blending surface. or gouging of the 
underlying surfaces where the linkage curves are closer together. The patches. linkage 
curves and directional curve are the same as those used in Figure 5.2. The blending 
surface is tangent plane continuous across both linkage curves. with q; the scalar bias 
functions as in equations (4.27) and (4.30), ; — 1.2 referring to the corresponding 
underlying surfaces and linkage curves. An alternative biasing method is a function 


based on the relative distance between corresponding points on the linkage curves. 


(LIEGE) - EXE. 


q; b. 9a) 


This formulation automatically accounts for changes in relative linkage curve position. 


vet allows for user control of the blend shape through variation of the factor. b,. The 
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Figure 5.4: Distance function bias factor, 6; = 2.5 


scalar factors b; are defined by the user for each linkage curve independently, and 
can themselves be constant values or functions of the parameterization of the linkage 
curves. An example of a blending surface using the distance formulation is shown 
in Figure 5.4. In Figure 5.3, the scalar bias function is q) = q2 = 0.4, with user- 
defined factors b; specified such that q; = constant. In the distance formulation of 
Figure 5.4, 6; = b, = 2.5 and the bias functions as given by equation (5.1) are now 
qı = q = 0.4|E,(£,) — E2(€;)|. The additional distance term adjusts the magnitude of 


q; to prevent the unwanted intersection of the underlying surface shown in Figure 5.3. 


Blending operations are local surface operations which must be carefully controlled 
to ensure the global geometry of the underlying surfaces is not changed. As noted in 
Chapter 1. the bias function gives the designer a measure of control over the interior 
of the blending surface. If the magnitudes of the bias factors are not appropriate. 
the blending surface can gouge one or both of the underlying surfaces. We give four 


examples of the effect of bias factor magnitudes. 


I 


O! 





The first surface is a one quarter of a cvlindrical surface expressed as a rational 


Bézier surface patch of orders 3x2. The knot vectors are 


25 0:0.0: 1; L1) 722020211) 


The control polygon vertices are 


Poo uL) ICD. 2-1) 
l — | 
Pio = (y 2 VOR: v Pii = E V2. ven ys! 


The linkage curve is a linear spline and has the knot vector and control points 
EN ee RE 05,10) P, = (0.5,0.0) 


expressed in the parameter space of the surface. The second surface is an integral. 


bilinear planar B-spline surface with knot vectors 
cun Ts peut) 
and control polvgon vertices 


Poo = DIE). 0) Po1 = D Pio = Kor) Pij — ‘Coe 2) 


The linkage curve is a linear spline expressed in the parameter space of the surface 


with knot vector and control polygon vertices 


ESOS UPS ORO PB, = (0.7, 0.0) 
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Figure 5.5: Curvature continuous blend. bias factor 6; = 6 


The directional curve is a linear spline with knot vector and control polygon vertices 


w 


V x 0. d O 


bo 


0,2) 1h = (2 


In the first example we show a well-specified blending surface. The bias factors 
are equal to 6; = 6 along each linkage curve with curvature continuity specified. The 
blending surface is shown in Figure 5.5. The bias factors are then increased to 6; = 12 
for each surface. Recalling the distance function as given in Equation 5.1, a larger 
bias factor weakens the biasing effect. We will call this a weak bias, as compared to 
a small bias factor which gives a strong bias condition. This is clearly demonstrated 
in Figure 5.6. For very strong bias factors, especially with curvature continuity 
specified along the linkage curve, it is visually difficult to tell the difference between 
the underlying surface and the blend. Furthermore. the stronger the biasing, the 
further the associated boundary effect is carried into the interior of the blend. If 


the bias factor is too strong, the blending surface will intersect the other underlying 
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Figure 5.6: Curvature continuous blend. bias factor b; = 12 


surface. In Figure 5.7, we try to pull the blending surface towards the planar surface 
by specifying an imbalance between the bias factors. The bias factor for the plane 
is by = 3, while the bias factor for the cylindrical surface is 6; = 10. The cylindrical 
surface is gouged by the blend because the biasing at the linkage curve on the plane 
is too strong. We found that to eliminate this tvpe of gouging, this biasing must be 
reduced. 

Increasing the biasing associated with the gouged surface just exacerbates the 
problem. creating not only inappropriate intersections of the blending surface with 


both underlying surfaces but self-intersections of the blending surface as well (see 


Figure 3.8). 


5.3 Marine Propeller Fillets 


Modern marine propellers represent particularly complex examples of sculptured sur- 


face engineering. Demanding performance criteria require accurate machining and 
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Figure 3.8: Blend with loop, bias factors 5; = 3 
19 





high surface quality. Hand finishing techniques are usually inadequate. As noted in 
Chapter 1. the use of numerically controlled machining requires explicit mathemat- 
ical representation of the entire surface, including the propeller blade-to-hub fillet. 
Failure to properly represent the fillet may result in removal of fillet material during 
blade machining, or gouging of the blade faces or hub may occur. 

The propeller fillet is both a structural and hydrodynamic surface addressing 
concerns which arise at the intersection of the propeller blade with the hub. The 
fillet contributes to the structural strength of the blade and aids in the reduction of 
vibrational blade loading [Allen 90]. By replacing a sharp edge intersection with a 
smooth transition surface. adverse hvdrodvnamic effects are also reduced. 

In Chapter 3. it was noted that parametric B-spline surfaces are ideallv suited 
for the representation of complex surfaces. Creation and interrogation of marine 
propeller surfaces represented with NURBS surfaces is presented in [Tuohy 91c] and 
¡Patrikalakis 90a]. A brief introduction to propeller geometry is provided in Ap- 
pendix A. 

A blending surface for the propeller blade-to-hub fillet will be constructed using 
two techniques. The first method uses the blending surface algorithm demonstrated 
earlier in this chapter. The second method is an automatic blending surface designed 
to replicate a standard propeller fillet geometry called the 3T-T/3 fillet [Gorman 79] 
¡Sabol 83]. A detailed discussion of the geometric requirements and implementation 
is included in Appendix B. The blending surface resulting from this implementation 
is a truly rational B-spline surface. For comparison, tangent plane and curvature 
continuous integral B-spline fillet blends are created using the general algorithm with 


linkage curves defined by the 3T-T/3 method. These blends are compared with the 


rational B-spline surface. 





5.3.1 Integral B-Spline Surface Fillet 


As the first example of the blending surface method applied to the propeller fillet 
problem, an integral B-spline blending surface is created between a B-spline surface 
representing a propeller blade and a rational B-spline surface representing a hub. A 
wire frame of the propeller blade. which has a maximum radius of ten. is shown in 
Figure 3.9. Thehub is a linear by quadratic NURBS surface representing a cylinder ot 
radius one. The linkage curve on the propeller is an isoparametric curve corresponding 
to v = 0.2. The v = 0 isoparameter curve of the propeller blade is selected as 
the directional curve. To find the linkage curve on the hub. the directional curve 
is projected orthogonally onto the hub surface. The linkage curve is the geodesic 
offset of the projection, with an offset distance of 0.2 units in three-dimensional 
space. The linkage curve is expressed in the parameter space of the hub. All of 
the curves were constructed interactively using Praziteles. Two blending surfaces 
were created using these surfaces and linkage curves. Figure 5.10 shows a tangent 
continuous blending surface with a bias factor b; — 2.0 applied along both linkage 
curves. lhe position tolerance specified for the blending surface approximation is 
0.01 The achieved tangent angle error from equation (4.51) is 0.2 degrees. Figure 5.11 
shows the same blending surface with propeller blade and hub. 

The blending surface in Figures 5.12 and 5.13 is curvature continuous with a bias 
factor b = 7.0 along the blade linkage curve, and b; = 9.0 along the hub. To provide 
an indication of the performance of the blending algorithm. a series of example blend- 
ing surfaces with these bias factors were created with different accuracy requirements. 
The results illustrating achieved position, tangent angle (equation 4.51), and relative 


curvature error (equation 4.56), are shown in Table 5.1. 


5.3.2 3T-T/3 Fillet 


The most common method for the design of a propeller blade-to-hub fillet is the 


3T-T/3 fillet developed at the Philadelphia Naval Shipyard. This design has been 
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Figure 5.9: B-spline representation of a propeller blade 








Figure 5.10: Tangent continuous blade-to-hub fillet blend for blade in Figure 5.9 
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Figure 5.11: Tangent continuous fillet between blade in Figure 5.9 and hub 





Table 5.1: Curvature continuous fillet with isoparametric linkage curve of Figure 5.13 





Required | Number o 
Position | Control pts Angle | Curvature 
Tolerance Position 
AAA 00009] oma] 630] 
oo] 3 [6 [0000| 0088] 61] 


[oo| s1 | 6 | 0.000130] 0.058] 651] 
[0.0001 [121 15 [00000422 0.081 6] 






























Required relative curvature tolerance=5.0% 
Required tangent angle tolerance=2 degrees 


Number of 
Position occa E. 
Tolerance S 
NE om ea 
| —9005| 5i | 8 | 0300196 | —0242| — 455] 
[— 0001 | 77 | 6 | 0300130 | — 9.58 | — 539] 


[0.0005 | 55 | e | 9.000130 | oos] 500] 
[0.001 [128 | 6 [0.000222] 0.0231 197] 







Required 




















Figure 5.12: Curvature continuous blade-to-hub fillet blend for blade in Figure 5.9 
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Figure 5.13: Curvature continuous fillet between blade in Figure 5.9 and hub 
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Figure 5.14: 3T-T/3 fillet blend for blade shown in Figure 5.9 


adopted by the Naval Sea Systems Command (NAVSEA) as standard geometry for 
naval propeller fillets [Sabol 83]. Experience has shown that the 3T-T/3 fillet is an 
effective transition both structurally and hydrodynamically. In Appendix B. the pro- 
cedure for constructing linkage curves and cross-link curves is presented in detail. 
The fillet surface is approximated as described in Section 4.4. The fillet is tangent 
continuous across the linkage curves. Figures 5.14 through 5.16 show the resulting 
3T-T/3 fillet surface. The blade in this example is the same as that used in Sec- 
tion 5.3.1. The hub is a linear by quadratic surface of radius 1.25 units. The shape 
function (Equation B.6) is specified to create a fillet with maximum thickness near 
mid chord, and 30% of the maximum thickness at the leading and trailing edges. The 


shape function implemented approximates a graph given in [Sabol 83]: 
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Figure 5.15: 3T-T/3 fillet blend for blade shown in Figure 5.9 
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Figure 5.16: 3T-T/3 fillet with blade and hub 





The mapping of the parameter u of the intersection curve of the blade to the hub to 


the parameter x in the shape function z - z(u) is 


"DS du (5:3) 
0505, => x= 4(0.5 — u) 53 
Mao de 05) (3.4) 
el > 7 = 4(1 — wu) 5 


The implementation of the 3T-T/3 algorithm is still experimental. As noted in 
Appendix B. the only error check currently implemented is position continuity ot 
the blending surface with the linkage curve on the hub. For the fillet illustrated in 
Figure 5.16. the achieved maximum position error is 0.03. 

To illustrate the robustness of the algorithm, a 3T-T/3 blend was created using 
a propeller blade of significantly different shape. The result is shown in Figures 5.17 


through 5.19. 


5.3.3 Pseudo 3T-T/3 Fillet 


The 3T-T/3 fillet was specifically designed to provide a suitable hydrodynamic and 
structural transition for marine propeller blades. We have no reason to believe the 
fillet created in Section 5.3.1 has similar properties. For the last two examples. we 
create propeller fillet blends using the method introduced in Chapter 4. with linkage 
curves of the 3T-T/3 blend. We will call this blending surface a pseudo 3T-T/3 fillet. 
In the first example. a tangent plane continuous fillet is created. This blend is most 
similar to the standard 3T-T/3 fillet since it. too. is tangent plane continuous. The 
fillet is illustrated in Figures 5.20 and 5.21, using the same propeller blade and hub 
as was used in Section 5.3.2. A series of examples were created with varying degrees 


of position accuracy specified. The results are given in Table 5.2. The bias factor b; 


for both the blade and the hub is 3.0. 
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Figure 5.18: 3T-T/3 fillet for blade in Figure 5.17 
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Figure 5.19: 3T-T/3 fillet surface for blade shown in Figure 5.17 


eure 5.9 


7 


.20: Tangent continuous pseudo 3T-T/3 fillet for blade in Fi 


gure 9 


r1 


30 





a, 
Mn RT 
AAN 
TOR 
Eg 
NS Dep 
| O 


V 
I» 
7 


id ES | 
na SS y ee 
CiT a AÁ— [ li S 
ca: SAO » | | 
448 LE a0 E ees N DM E | 
MA SS ARICA Neu I 
A ASS INS IN 
AES EA II INN LIN 
pilo e o EE CROYDON | 
SES SSI LN | 
d RI, IN N ST N SAILS SSE NEAT XT 
/ REESE SSN 
ERS SSS SNS} 
N SINN NN | NS N 8 | SAR 
PSS SSA DIN 


I 
T 
Af 

CLL LL 


v 


) 


| Ln 
Ly 
TH TA 


[m 
I d 


NI) 


4 


pi 
N 
| 


7 
Ó 

| W 
ey 
A 
7 


TY 


P 
| 
CIS 


/ 
f 
/ 


Ew 


j 
! 
OSA y N 


E 


E H 
If} 
TEE 
LETT | 
ij IT 177 77 
MWA 
Uil 
ENS - 
SY 
SS 


/ 


Figure 5.21: Tangent continuous pseudo 3T-T/3 fillet for blade in Figure 5.9 


Table 5.2: Tangent continuous pseudo 3T-T/3 fillet 


| Required tangent angle tolerance=2 degrees | 


Required | Number of 

Tolerance u | v | Position TÄnglerdes)] 
[oo | 3 | 90s] 135] 
ns | 1 Foo — 033. 
0.001] 65 1 | (30860 | — 0:5] 
[90005 [ 5 | 1 | 0300113 | or] 
[  6900L [109 | 1 [9.000482 | — X067] 
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Table 5.3: Curvature continuous pseudo 3T-T/3 fillet 





Required relative curvature tolerance=10.0% 
Required tangent angle tolerance=2 degrees 
















Required | Number of Achieved Errors 


Position pem Le Angle NN 
Position | (degrees) | (Percent) 








Tolerance | u | 


See n 
0.051 50 [6 | 0.0ossı | 0.517 — 939] 
0.011 65 | 6 | 0.000161 | 0200] 6.10] 
00005 | 5 | 6 | 0000114] 0266] — i12] 
[— 6300. [109 | 6 |0.0000462| —.06:0| — 5i 
| .00005]128 | 6 [0.000188 | 0.0584] 201] 






Required relative curvature tolerance- 1.0706 
Required tangent angle tolerance=2 degrees 


Number of Achieved Errors 
Position 


P rol m Angle | Curv um 
Tolerance | u | Position | (degrees) | (Percent) 


TA 
0051 71 [5 | 0.000556 | 0.0644] — 0.901 | 
| 000| 19 | 6 | 0.000011 | —.058L| — 0.950 | 
[ 0.0005 | $6 | 6 | O.000IÍ| 0.0825 | — 0.956 | 
[0.0001 [116 | 6 [0.000462 | 0387 | — 0.165] 
| 9.00005 [130 | 6 [0.000088 | 0.0538] 0.801 | 


















Required 






To complete the examples, a curvature continuous pseudo 3T-T/3 blend is illus- 
trated in Figures 5.22 and 5.23. The bias factor for the blade is 6, = 7.0, and b» = 9.0 
tor the hub. These bias factors were selected experimentally to create a blend which 
most resembles the standard 3T-T/3 fillet. As with the other examples. a series of 
blending surfaces was created with varying degrees of accuracy. The results are shown 
Bie lable 5.3. 

Review of the tables indicates the number of control points required to represent 


a curvature continuous blend is heavily dependent on the specified position tolerance. 








For a required position tolerance ot 0.01. 10 additional control points in the u direction 
are required to create a blend with a relative curvature error tolerance of 1096. when 
compared with the tangent continuous blend listed in Table 5.2. Decreasing the 
relative curvature error tolerance by one order of magnitude for the same position 
tolerance of 0.01 more than doubles the required number of control points. If the 
required position tolerance is now reduced to 0.001, the curvature continuous blend 
with a relative curvature error tolerance of 10% has exactly the same number of 
control points as the tangent continuous blend. Decreasing the relative curvature 


error tolerance to 1% requires only 14 additional control points. 
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Chapter 6 


Conclusions and 


Recommendations 


In this thesis. a general algorithm for the representation and approximation of 
blending surfaces between two non-uniform rational B-spline (NURBS) surfaces using 
a surface of the same class is presented. As discussed in Chapter 1, blending surfaces 
are generally procedural surfaces which are difficult to process. Further, there are no 
standard data exchange formats for procedural surfaces. Our use of NURBS surfaces 
is Important because it facilitates the exchange of geometric information between 
various geometric modelling systems, and design and manufacturing svstems. 

Ihe method approximates the blending surface with an integral B-spline surface 
which is up to curvature continuous with the underlying surfaces along the linkage 
curves. Several methods are available for specifying the linkage curves. The method 
allows the degree of continuity along the linkage curves to be specified independently. 
enhancing its flexibility. Bias functions are control parameters which influence the 
shape of the blending surface. We demonstrated that if weak biasing is used. the ef- 
fect of the boundary condition is confined to a narrow region along the linkage curve. 
Strong biasing propagates the boundary conditions far into the blending surface in- 


terior. If the biasing is not carefully controlled. however. inappropriate intersections 
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with the underlying surfaces can result. Our implementation uses a bias function 
based on the relative distance between corresponding points on the linkage curves. 
This technique adjusts the magnitude of the bias so that in regions where the linkage 
are close together, the probability of gouging the underlying surfaces is reduced but 
not eliminated. Careful user control is required. 

An additional control element is the directional curve. Our implementation guar- 
antees the cross-link curve will be oriented such that its tangent direction at the 
endpoints is towards the directional curve through the use of an angle comparison 
algorithm. 

The linkage curves are approximated with low degree integral B-spline curves. 
An adaptive sampling technique adds knots at appropriate locations until the ap- 
proximation meets all specified tolerances. The cross-link curves are also represented 
with low degree B-spline curves. The cross-link curves are constructed such that the 
surface will satisfy the specified continuity conditions at the linkage curves. Meth- 
ods are implemented for position, tangent plane and curvature continuity along the 
linkage curves. The cross-link curves then form a family of curves which become the 
isoparameter curves of a lofted blending surface. The procedure again uses adaptive 
sampling to improve the approximation. The surface is sampled at the midpoints of 
the common knot vector formed by merging the knot vectors of the linkage curves. 
The error checking routine compares the position and surface normal vector of the 
blending surface with the same information from the corresponding underlying sur- 
face. To compare curvature, normal curvature at the endpoints of an isoparameter 
curve extracted from the blending surface is compared with the normal curvature 
of a cross-link curve created at the same node using the underlying surfaces in the 
normal manner. Points which fail to meet the prescribed tolerances are added to the 
knot vector and the surface relofted. Throughout the improvement process. position 


accuracy is maintained. 


We illustrate the practicality of the method by creating marine propeller blade-to- 
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hub fillets. Three types of fillets were created with varying degrees of continuity, and 
each was evaluated over a range of specified position accuracies and. where applica- 
ble. relative curvature tolerances. Numerous examples of fillet surfaces with varving 
degrees of accuracy were examined. We found that required position tolerance with 
curvature continuity specified has a significant effect on the resulting surface. When 
the position tolerance is relatively course (of order 0.01), the magnitude of the rela- 
tive curvature tolerance controls the number of control points necessary to represent 
the fillet. In the case of the pseudo 3T-T/3 fillet. decreasing the relative curvature 
tolerance by one order of magnitude nearly doubles the number of required control 
points. If the position tolerance 1s smaller. the change in curvature tolerance has 
less effect. It was demonstrated that position accuracy is maintained during the loft- 
ing process, with the achieved position error one to two orders of magnitude better 
than the position tolerance specified. It was also noted that there is little difference 
in terms of data storage between representations of tangent continuous fillets and 
curvature continuous fillets when created with small position tolerances. 

The shape of the linkage curves and the underlying surfaces has an impact on the 
achievable degree of continuity. We found that the fillet designed with an isopara- 
metric linkage curve on the blade could not achieve better than 4% relative curvature 
error, no matter how small we made the required position tolerance. In comparison. 
it was possible to create fillets with less than one percent curvature error when using 
the 3T-T/3 linkage curves. even with a relatively coarse position tolerance of 0.01. 

Ihe pseudo 3T-T/3 fillet effectively illustrates the practical application of the 
algorithm. The rational B-spline 3T-T/3 fillet has the advantage of being automatic 
in the creation of the fillet. On the other hand. the general algorithm can create 
a fillet which approximates the rational fillet while retaining the flexibility of user 
control through the bias functions. Furthermore. the fillet can be either tangent 


plane or curvature continuous while the rational fillet. by its design. can be at best 


only tangent plane continuous. 
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The blending algorithm presented in this thesis has the following advantages: 


l. The approximation of the linkage curves and blending surface is adaptive. Sam- 


pling and knot addition depends directly on the distribution of error. 


2. The blending surface approximation is in a standard format allowing direct 
exchange between different geometric modelling systems which use the standard 
data exchange formats [Smith 88]. Further, the surface is of low degree which 


permits processing on less capable geometric modelling systems. 


J. Ihe user has direct control of the shape of the blending surface through speci- 


fication of boundarv conditions. bias functions. and the directional curve. 


The following areas are identified for further work: 


l. Extend the implementation of the method in Chapter 5 bevond cross-link curve 


curvature continuity to cover normal section curvature continuity. 


2. Provide interrogation methods for the blending surface to ensure there are 
no unwanted intersections with the underlying surfaces or self-intersections 


[Kriezis 90]. This might be included as part of the adaptive sampling process. 
3. Generalize the algorithm to allow the blending of more than two surfaces. 


4. Experimentation with the directional curve to gain greater control of the blend 


shape. Experimentation with alternate parameterization schemes of the linkage 


curves and the resulting blending surface. 


5. Improving automation of methods for controlling shape and position of the 
linkage curves. This might include improving the accuracy and efficiency of the 


rational 3T-T/3 algorithm currently implemented. 


6. Automation of bias function specification to eliminate unwanted intersections. 





T. Improve methods for visualizing continuity conditions across the boundaries of 


blending and underlying surfaces. 
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Appendix A 


Marine Propeller Geometry 


Marine propeller blade geometry is usually described with reference to a Cartesian 
coordinate system. The orientation presented here is as suggested by [Kerwin 86]. 
The z-axis is aligned along the axis of rotation of the propeller and assumed to be 
positive in the direction of mean flow. The y-axis is vertical and passes through the 
mid-chord of the root section of one blade. The z-axis is oriented to form a standard 
right-handed Cartesian coordinate system (see Figure A.1). A cylindrical coordinate 
system with r the radial coordinate from the z-axis and @ the rotational coordinate, 
measured about the z-axis such that 6 = 0 on the y-axis. is also used. 

The propeller blade is typically constructed of a number of two-dimensional foil 
sections lying on the expanded planes of cylindrical surfaces coaxial with the z-axis. 


The chord of the section is the two-dimensional straight line segment connecting the 





leading and trailing edges. In three dimensions the chord is an arc along a helix and 
is a measure of the pitch angle of the blade at the particular radial distance of the 
section. The leading edge is determined by advancing, in the negative direction, along 
a cylindrical or conical helix defined by the radially distributed pitch angle a distance 
equal to half the chord length. The shape of the foil section is defined in terms of 


the chordwise distribution of camber, f(s), and thickness, t(s) (see Figure A.2). The 





parameter s is a coordinate defined along the chord. generally 0 at the leading edge 
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Figure A.1: Propeller blade geometry (from [Kerwin 86)) 
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EDGE 


Figure A.2: Foil section 
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and | at the trailing edge. The line traced by f(s) is the camber line, defined as the 
locus of the midpoints of straight line segments spanned across the section such that 
the normal to the segment at its midpoint is tangent to the camber line. [The camber 
function defines the shape of the camber line relative to the chord of the section. The 
half-length of the segments which span the section is defined by the thickness function. 
These segments define pressure and suction faces of the blade. The pressure face is 
the blade surface facing aft which experiences an increase in pressure when propelling 
a ship ahead. while the suction face, or back, experiences a decrease in pressure. The 
pressure differential results in lift. The root section is the foil section formed bv the 


intersection of the blade with the hub. prior to any blending or filleting operations. 
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Appendix B 


Creation of a 3T-T/3 Marine 
Propeller Fillet 


This Appendix presents a detailed description of our implementation of the stan- 
dard 3T-T/3 marine propeller fillet using the general algorithm of this thesis. 

First, we recall the general method for the creation of a blending surface. We 
specify linkage curves on the blade surface and the hub. In place of a directional curve 
we construct the cross-link curves in a fillet plane defined by a geometric function. 
Ihe cross-link curves are a family of rational B-spline curves which satisfy another 
specific geometric relationship. Finally, the fillet surface is created by lofting the 
cross-link curves. 

Each cross-link curve is a planar curve constructed of two circular arcs. One of 
the arcs is tangent to the propeller blade, and one is tangent to the hub surface. In 
addition, the arcs are tangent to each other. The circular arc tangent to the blade 
has a radius of 3T, while the arc tangent to the hub has a radius of T/3. hence the 
name given the fillet. The dimension, T', is defined in section B.1. below. Figure B.1 


illustrates the fillet geometry. The generation of the fillet is as follows: 


e Compute the intersection curve of the blade and hub. and the maximum thick- 


ness of the root section. 
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Figure B.1: 3T-T/3 fillet cross-link curve geometry 
Find the linkage curve on the blade surface. 


Establish the fillet plane orientation and the center of the circular arc associated 


with the linkage curve on the blade. 
Compute the intersection of the fillet plane with the hub. 


Find the center of the T/3 circular arc and its associated contact point on the 


hub. This is the linkage point on the hub. 


Construct circles in the fillet plane using the computed centers. Subdivide the 


circles at the points of tangency with the underlying surfaces and with each 


other. 

Join the arcs to form the cross-link curve. 

Approximate the blending surface from the family of cross-link curves by lofting. 
Sample and improve the approximation. 
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B.1 Linkage Curve on the Blade 


The linkage curve on the blade is defined to lie on the blade surface at an arc 
length distance T' from the blade-to-hub intersection. The original development of the 
3T-T/3 fillet suggests that the dimension T' is the helical arc length in the root section 
from the pressure face to the suction face of the unfilleted blade contours. Ihe di- 
mension is a local arc length parameter in that it varies with chord length [Sabol 83]. 
NAVSEA further recognizes the imprecision of the definition and so provides manu- 
facturing guidance to modify the fillet gradually around the section so that it is fair. 
[Sabol 83] suggests that T can be defined as a function of the maximum thickness of 
the root section and a shape factor which varies with chord length. 

The maximum thickness of the root section of the blade is found using techniques 
for feature extraction from B-spline representations of marine propellers presented 


by [Patrikalakis 90a]. The propeller blade is represented as a bicubic non-uniform 


integral B-spline surface: 


n=1m-1 


R. (uv) = > 3 peu sU) (B.1) 
i=0 j=0 

The parameter u corresponds to the chordwise direction from the trailing edge on the 
suction face to the leading edge, to the trailing edge on the pressure face. The v pa- 
rameter corresponds to the spanwise direction, hub to blade tip. The v-isoparameter 
lines thus approximate the intersections of the blade surface with cylindrical surfaces 
coaxial with the propeller axis of rotation. Without loss of generality. we assume a 
cylindrical hub. and the section of interest 1s the root section formed bv the intersec- 
tion of the cylindrical hub with the blade. Since the intersection curve is known to be 
a simple regular curve without self-intersections. we use a lattice evaluation technique 
to find the intersection. Data points representing the intersection of the cylindrical 
hub surface with isoparameter curves in the v direction (i.e.. lines of constant u). 


extracted from the blade surface hub to tip, are calculated using a simple bisection 
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method. While more elaborate or efficient methods are available. this one is quaran- 
teed to find the appropriate intersection. The intersection curve i(t) is approximated 


as a cubic integral B-spline curve 


n=l 
iS x I; B, 4(t) (B.2) 
1=0 
with knot vector T = (tg,ti,..-,tn4n—-1), and N = 4. The approximation of the 


* od 


intersection data is similar to the method described in Section 4.2. In this case. (B.2) 


interpolates the intersection data at the nodes £j; of the knot vector 


1 3 
i= a ditiek 0Sisn-1 (B.3) 
> k=] 
such that 
n—1 
ee ec) =e), 0S<:<n-1 (B.4) 
9-0 


Since (B.4) is of the same form as (4.4) it can be written in matrix form. 


(NJ = {i}, (B.5) 


and, like (4.5), can be solved using standard linear algebra techniques. The approx- 
ımatıng B-spline ıs expressed ın terms of normalized cylindrical coordinates. The 
approximating intersection curve is virtually expanded to a plane by cutting it along 
a generatrix not intersecting the section. The section curve is thus a planar curve 
assumed to be at least C? continuous. lying on the expanded cylinder. To find the 
camber line, a correspondence is established between parameters on the pressure and 
suction faces in accordance with the definition of camber line as given in Appendix A. 
These are referred to as matching parameters. [Patrikalakis 90a] develops a system 
of differential equations which. when integrated with respect to the arc length of the 
section chord. completely define the camber line. The integration proceeds until an 


irregular point is reached. An irregular point represents a point where the tangent 
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vectors to matching parameters are parallel to each other. If the segment connecting 
the matching parameters is also normal to the tangent vectors, the point is called 
the critical point. |Patrikalakis 90a] shows that there is only one critical point on the 
camber line and that 1t corresponds to the maximum thickness of the section. Thus. 


to find the maximum thickness of the section we use the following procedure: 


1. The root section curve is found and approximated as an integral B-spline curve. 


I(t). 


2. The system of differential equations given in [Patrikalakis 90a] is integrated us- 


ing a variable step Adams method [NAG 89] until an irregular point is reached. 


3. A Newton-Raphson method is then used to accurately solve a nonlinear system 
of equations to find the critical point. The maximum thickness. T5,4, of the 


section is derived directly from the definition of the critical point. knowing the 


section curve. 


Readers interested in the details of the development of the systems of equations are 
referred to [Patrikalakis 90a]. 

The parameter 7’ can now be found. As noted earlier. T' is a function of the maxi- 
mum section thickness and a shape factor. The shape factor is defined to produce the 
maximum fillet thickness near midchord and a minimum thickness near the leading 
and trailing edges of approximately 1/3 the maximum thickness. In our implementa- 


tion the shape factor is defined as a general cubic function of the form 
f(z) = Ar? + Br? +Cr+D=0. (B.6) 


where 0 € x € 1. We define a function z — z(u) (equations (5.2) to (5.5)), to establish 
a mapping between the variable z and the parameter u of the propeller blade-to-hub 
intersection curve, expressed in the parameter space of the blade. to achieve the proper 


shape of the fillet. The local thickness dimension along the intersection curve is then 
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defined as 
Ti(u) = Tnas f(z(u)) (B.7) 


To construct the linkage curve on the blade, we evaluate the approximating 
B-spline curve of the intersection curve. At each value of u. we extract an isopa- 
rameter curve, u = constant. The corresponding point on the linkage curve is found 
along this isoparameter curve, measured in terms of arc length the distance Tı(u). 
The linkage curve I’,,op 1s approximated in the parameter space of the propeller blade 
as before, using as the initial knot vector the knot vector of the approximation to the 
intersection curve. 

The general algorithm requires two linkage curves with each curve defined inde- 
pendently. In the case of the 3T-T/3 blend, however, there is a specific geometric 
relationship between a point on the linkage curve of the blade and the associated 
point on the linkage curve of the hub. Because of this dependency, we do not specify 


the linkage curve on the hub directly but find the corresponding points on the hub 


while constructing the cross-link curves. 


B.2 Fillet Plane 


To establish the fillet plane, we evaluate the approximations to the linkage curve on 
the blade and the blade-to-hub intersection curve at the nodes £; of the common 


knot vector. At each node there is a unique fillet plane described by the equation 


[Hildebrand 76] 
a Dy FCD =O: (B.8) 


The fillet plane is established by the normal vector to the blade surface computed on 


the intersection curve at a node. and the corresponding point on the linkage curve 


(see Figure B.2). 


Ihe unit normal vector to the blade surface, Nyrop, at a node of the intersection 
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Figure B.2: Fillet plane orientation 
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curve, is found with equations (4.22) to (4.24). A second vector p is defined as 


) 


crece 
[Eoo (£i) y» I(£&;)| 


pu (B.9) 


where E,,,,(6;) 1s the linkage curve position on the blade evaluated at the node £;. 
and I(£;) the corresponding point on the approximation to the intersection curve. Let 


n, be a unit vector which is normal to the fillet plane 


n (€;) = p(&) ^ He (B.10) 
Á Ip(£;) ^ nic 
In terms of equation (B.8). 
n, = Ai 4- Bj - Ck. (B.11) 
The remaining term in (B.8) is 
D = prop se) Í sale.) (B12) 


B.3  Cross-Link Curves 


To construct the cross-link curves, the centers of the two circular arcs shown in 


Figure (B.3) and the three points of tangency must be determined. 


B.3.1 Centers of Circular Arcs 


Ihe center of the 3T circle is found first. The first parametric derivative. R* = 
R3 4(u.v), of the propeller blade in the v direction (3.21) is computed at the corre- 


sponding point on the linkage curve, E,,,,(£;). The vector q defined by 


GE mn, AR. (B.13) 
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Figure B.3: 3T-T/3 detailed geometry 
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gives the direction to the center of the 3T circle. L. 


Mi +3 (B.14) 


where Q(u,v) is the point on the blade surface corresponding to the point Eyrop(£i) 
on the linkage curve. expressed in three-dimensional coordinates. Likewise. the vector 
q, defined in Equation (B.13), is evaluated at the corresponding point on the linkage 
curve. 
To locate the center of the T/3 circle. S. we must find the intersection of the fillet 
plane with the hub. The hub is given as a rational B-spline surface. 
Eoo usu Pi BS ()B5s(v) 


Hxy3(u.v) = een B.15 
EUN DE x UI a) P 


Mere H(u,v) 2 (H,, H,, H,), and 


Sy a E v) 


(có 2.5-o Ui Bi v (u) Bja(v) 


H,(u.v) = (B.16) 


H, and H, are defined similarly. The intersection of the fillet plane defined by (B.8) 


is therefore 
eo deer aD = 0, (o d) 
Since A. B, and C are constants from (B.11), (B.17) can be rewritten as 
n—1 8 
Due v) = D > Eire + JB + Cee — D)Bi x(u)B;a(v) = 0 (B.18) 


¿¡=0 ¡=0 


EE = O, then df = f,du + f,dv = 0, or 


du f 


PE a (B.19) 


This formulation requires special treatment when f, = 0 or f, = 0. The difficulty can 


be overcome by expressing the intersection curve as a function of the parameter. t, 


which leads to the following system of differential equations [Kriezis 90] 


df du dv 
eee ace OO _ B.2 
n if pues cd (B.20) 
3i EROS 
DENT Ju (B.22) 


4 7 men 


The system of differential equations (B.21) and (B.22) are solved using an Ádams- 
Bashforth technique [NAG 89] to find points along the intersection curve of the fillet 


plane with the hub surface. The geometric requirements for the T/3 circle are that it 


l. be tangent to hub, and 


2. be tangent to the 3T circle at a point between the points of tangency of the 


circles on the hub and blade. 


Using simple geometric construction it can be shown that the center of the T/3 circle 
is located at the intersection of a circle of radius (3T; — Tj/3), centered at L. and a 
curve 1n the fillet plane offset from the intersection curve of the fillet plane with the 
hub at an offset distance 7;/3 (see Figure B.4). The center of the circle must be in 
the fillet plane and, therefore, satisfy (B.8). The second requirement is that it lie on 


a circle of radius T;/3. centered at L. The equation of a sphere. centered at L with 


radius (3T; — T}/3) is 
2 NA E / I 2 o: 
EL + == 0) +(2-—() = i (B 23) 


where (€,4,C() are the three dimensional coordinates of L. The intersection of the 
fillet plane with this sphere is a circle of radius (T, — T,/3). We can therefore find the 
center S by restricting the solution to points which satisfy the equation of the fillet 


plane (B.8). By using the equation for a sphere we free ourselves from working with a 
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Figure B.4: Center of T/3 circle 


local coordinate system. We compute the tangent vector t(t) to the fillet plane to hub 
intersection curve by solving the system of differential equations (B.21) and (B.22). 
The vector n* is the vector product of the normal vector to the fillet plane and the 


tangent vector, t 


ur mtt. (B.24) 
Ihe center of the 7;/3 circle is found by computing 


S- ins + Hm, (u(t), v(t)) (B.25) 


using an adaptive step size marching method. The solution of (B.25) which solves 


both equations (B.8) and (B.23) is the center, S. 


B.3.2 Generation of Cross-Link Curves 


We now have all the information necessary to create the cross-link curves. namely, 
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1. the center of the 3T circle. L = (Lz, L,, L,), 


2. the point of tangency of the 3T circle with the propeller blade. Q(€;) = 
(Qe, ale 


Bethe center of the T/3 circle, S = (Sy, S,, S2), 


Ir YyyJz 
1. the point of tangency of the T'/3 circle with the hub, and 
5. the normal vector to the fillet plane. n,. 


Each circle is defined as a rational B-spline curve 


8 DD 
A Zizo WiP Biat) (B.26) 
io w¡B;3(t) 


| Jm h;K;B;s(s) 
ara ee oe B97 
(s) BS eae 


where 0 < t,s < 1. Let A be the circle of radius 377, in the fillet plane with its center 
at L. and C the circle of radius 7)/3 centered at S. also in the fillet plane. We seek 
to subdivide these curves at the points of tangency with the underlying surfaces and 
with each other, then join the segments to form the cross-link curve. 

Let the point QO = (Qz, Qu, Qz) (equation (B.14)) be the coordinates of the point 
of tangency between A(t) and the blade surface. We define the parameters t,.s, 
as the values associated with the points of tangency, J, where the circles are to be 
subdivided. Given A(t) = (z(t), y(t), =(t)), the point of tangency of A(t) with the 


blade surface is the value tı such that 


(x(t1) — Qz)* + (y(t1) — Qy)? + (2(t1) — Q.)? = 0 (B.28) 


is satisfied. A minimization technique [NAG 89] is used to find t,. A similar expression 


is evaluated to find the parameter s, where C is tangent to the hub. 
We find the parameters associated with the intersection of the two circles in the 


same manner. In this case, we evaluate the intersection of each circle with a sphere 
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located at the other center. The circle described by A(t) lies on the sphere 
(x(t) — Lz)” + (y(t) — Ly)” + (=(t) — Lo) - BT) =0 (B.29) 


Similarly, the smaller circle C hes on the sphere 
(z(s) — $2)? + (y(s) — Sy)? + (2(s) — 5.)? — ($4)? =0 (B.30) 


At the point of tangency between the two curves. the spheres (B.29) and (B.30) are 
also tangent. Using a minimization technique with respect to t. we find the intersection 
of A(t) with the sphere given by equation (B.30). Similarly. the parameter associated 
with C(s) 1s found minimizing equation (B.29) with respect to s. 

Knots are added to the knot vectors of A (t), C(s), at the parameter values t;, s,, ) = 
1,2, and the curves subdivided using the Oslo algorithm described in Section 3.2.1. 
The circular arcs are rational B-spline curves of order 3, geometrically identical to 
the segments from the original circles between the parameters t,,t2 and sı,$s2, each 
with a new knot vector normalized over the interval [0.1]. Referring to the segments 
as A(t) and C(s) as before, the curve A(t) has m control points, while C(s) has n 
control points. The segments are now joined to create the cross-link curve. 

To join the two curve segments we combine the knot vectors and the control points. 
This results in a double control point at the intersection of the curve segments, one 
of which must be removed. Recalling the formulation of a rational B-spline curve, 

A Lizo wiPiBis(t) (B31) 
Lico WiBialt) 


n- LR. 
er Liso UK iBiats) (B.32) 
i=0 h;B;3(s) 


the double control point is Pp and K,_,. Because the homogeneous coordinates are 
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not. in general. equal. we normalize both curves: 


Eu SUELE) 
ENIM 0 ———— SS — i255 
) D w;hn-ıBislt) 
wi h;wok; B;3(t) 
=p Sate eee A B 4 
Cts) a Dal) aa 


where wo and h„_ı are the homogeneous coordinates of the double control point. We 
are now free to remove one of the control points at the intersection. The knot vectors 


are combined as follows. The knot vectors are 


D Een) (B.35) 


l 


(So, S1) +---Sn+2) (B.36) 


where J is associated with A, and S with C. To combine the knot vectors. we remove 
the multiple knots 54, 5441, 5442. Because we remove one control point. the knot to is 


also removed. Finally, the range of the knot vector T' is changed 


T (LOI O te FI) (B.37) 


to keep the knot vectors T and 5 from overlapping. The combined knot vector spans 


the interval [0, 2]. and is then normalized over the interval [0. 1]. 


B.4 Approximation of 3T-T/3 Fillet Surface 


The approximation of the fillet surface uses the same method as the approximation 
of the integral B-spline surfaces described in Chapter 4. The knot vectors of the 
individual cross-link curves are merged to form a common knot vector. The cross-link 
curves are then lofted to create a rational B-spline surface. 

To improve the approximation. we test for position continuity of the blending 


surface along the boundary of the fillet on the hub. We first find the intersection of the 
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T/3 circle with the hub as shown in Section B.3.1. for a cross-link curve constructed 


at the midspan between the nodes of the common vector. 


we 
t = e (B.38) 
The surface approximation error 1s 
EE EE v) 0 sm 1. Broo 


where H(t;. v) 1s the intersection point of the T/3 circle on the hub expressed in three 
dimensional space. and S(t,, v) is the corresponding point on the blending surface. If 
the error function 1s greater than the specified tolerance. e. the point t; is added to 
the knot vector. The blending surface is relofted with the enriched knot vector until 


the blending surface passes the error test at all points t,. 
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